AD 

TECHNICAL  REPORT  ARBRL-TR-02367 

RECONSTRUCTION  OF  A  BLAST  FIELD  FROM 
PRESSURE  HISTORY  OBSERVATIONS 

Aivars  Celmips 

September  1981 


US  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  COMMAND 
BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND,  MARYLAND 


! 


Approved  for  public  release;  distribution  unlimited. 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Secondary  distribution  of  this  report  by  originating 
or  sponsoring  activity  is  prohibited. 

Additional  copies  of  this  report  may  be  obtained 
from  the  Defense  Technical  Information  Center,  Cameron 
Station,  Alexandria,  Virginia  22314. 


The  findings  in  this  report  are  not  to  be  construed  as 
an  official  Department  of  the  Army  position,  unless 
so  designated  by  other  authorised  documents. 


Tha  use  jf  :raie  nones  or  menu, ft asrurera '  nones  in  this  report 
Joes  not  constitute  indorsement  of  any  commercial  product . 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER 

TECHNICAL  REPORT  ARBRL-TR-02367 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (and  Subtltls) 

Reconstruction  of  a  Blast  Field  from 

Pressure  History  Observations 

S.  TYPE  OF  REPORT  &  PERIOD  COVERED 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORf  a) 

Aivars  Celmi^is 

8.  CONTRACT  OR  GRANT  NUMBERfs) 

9  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Ballistic  Research  Laboratory 

ATTN :  DRDAR-BLI 

Aberdeen  Proving  Ground,  MD  21005 

10.  PROGRAM  ELEMENT,  PROJECT,  TASK 

AREA  &  WORK  UNIT  NUMBERS 

1L161102AH43 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Armament  Research  &  Development  Command 

US  Army  Ballistic  Research  Laboratory 

ATTN :  DRDAR-BL 

Aberdeen  Proving  Ground,  MD  21005 

12.  REPORT  DATE 

SEPTEMBER  1981 

13.  NUMBER  OF  PAGES 

89 

14.  MONITORING  AGENCY  NAME  &  ADDRESS  (If  different  from  Controlling  Office) 

IS.  SECURITY  CLASS,  (of  thla  report) 

UNCLASSIFIED 

IS©.  DECLASSI  F|  CATION/ DOWN  GRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  this  Report) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  ebetract  entered  in  Block  20,  If  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  reverse  aide  If  neceaaary  and  Identify  by  block  number) 

Blast  Field  Model  Fitting 

Overpressure  History  Data  Numerical  Integration 

Spherical  Blast  in  Ideal  Gas  Variance  Estimates 


20.  ABSTRACT  (Continue  am  reverse  atria  if  neceaaary  and  Identify  by  block  number)  dll 

This  report  describes  a  method  for  the  calculation  of  parts  of  a  spherical 
blast  field  from  pressure  history  observations  at  stations  placed  at  various 
distances  from  the  blast  center.  The  method  consists  of  a  determination  of  the 
pressure  field  by  model  fitting  and  of  a  subsequent  numerical  integration  of 
the  flow  governing  equations.  The  results  provide  a  complete  flow  description 
(pressure,  density,  and  particle  velocity)  as  well  as  estimates  of  the  accuracy 
of  the  results  within  the  region  of  observations.  A  check  of  some  basic 

(Continued) 

DD  1  JAM ^73  1473  EDf-nOM  OF  »  NOV  65  )S  OBSOLETE  UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Deta  Entered) 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGEfWhan  Data  Entered) 


20.  Abstract  (continued) 

assumptions  (spherical  symmetry  and  negligible  viscosity)  is  provided  by 
an  independent  calculation  of  the  particle  velocity.  Application  examples 
are  shown  for  a  theoretically  computed  strong  blast  field  and  for  data 
from  a  high  energy  blast  experiment. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfWhen  Data  Entered) 


TABLE  OF  CONTENTS 


Page 

LIST  OF  TABLES .  4 

LIST  OF  FIGURES .  5 

1 .  INTRODUCTION .  7 

2.  BASIC  ASSUMPTIONS  AND  THEORY .  8 

3.  NUMERICAL  INTEGRATION  AND  ACCURACY  ESTIMATES . 12 

4.  DETERMINATION  OF  SHOCK  FUNCTIONS  BY  MODEL  FITTING  .  18 

5.  ANALYSIS  OF  PRESSURE  HISTORIES . 23 

6.  DETERMINATION  OF  PRESSURE  FIELD  FUNCTIONS  BY  MODEL  FITTING . 26 

7.  A  THEORETICAL  TEST  EXAMPLE . 28 

8.  APPLICATION  TO  DATA  FROM  MISERS  BLUFF  II . 49 

9.  CONCLUSIONS . 80 

ACKNOWLEDGMENT . 82 

REFERENCES . . 

LIST  OF  SYMBOLS . 84 

DISTRIBUTION  LIST .  oc 


3 


LIST  OF  TABLES 


Table  Page 

1  Parameters  of  the  Theoretical  Test  Example .  29 

2  Assumed  Accuracy  of  the  Theoretical  Data .  29 

3  Shock  Parameters  of  the  Theoretical  Example  .  35 

4  Overpressure  Field  Parameters  of  the  Theoretical  Test  Example  .  39 

5  Ambient  and  Shock  Data  of  MISERS  BLUFF  II .  50 

6  Shock  Parameters  of  MISERS  BLUFF  II  .  51 


4 


LIST  OF  FIGURES 


Figure  Page 

1  Integration  Paths . ^3 

2  Computation  of  Flow  History  at  a  Given  Distance . 14 

3  Location  of  History  Observations  in  Theoretical  Test  Example  .  .  30 

4  Shock  Overpressure  vs.  Distance  in  Theoretical  Test  Example.  .  .  32 

5  Shock  Overpressure  vs.  Time  in  Theoretical  Test  Example . 33 

6  Shock  Distance  vs.  Time  in  Theoretical  Test  Example . 34 

7  Fitting  of  Pressure  at  120  m  in  Theoretical  Test  Example  ....  35 

8  Fitting  of  Pressure  at  170  m  in  Theoretical  Test  Example  ....  3-7 

9  Log-log  Plots  of  History  Parameters  vs.  Distance  in  Theoretical 

Test  Example . ^g 

10  Pressure  History  at  150  m  in  Theoretical  Test  Example . ^ 

11  Particle  Velocity  History  at  150  m  in  Theoretical  Test  Example  .  42 

12  Density  History  at  150  m  in  Theoretical  Test  Example . 43 

13  Dynamic  Pressure  History  at  150  m  in  Theoretical  Test  Example.  .  44 

14  Overpressure  History  at  160  m  in  Theoretical  Test  Example.  ...  45 

15  Particle  Velocity  History  at  160  m  in  Theoretical  Test  Example  .  46 

16  Density  History  at  160  m  in  Theoretical  Test  Example . 47 

17  Dynamic  Pressure  History  at  160  m  in  Theoretical  Test  Example.  .  48 

18  Shock  Overpressure  Versus  Distance  in  MISERS  BLUFF  II .  52 

19  Shock  Overpressure  Versus  Time  in  MISERS  BLUFF  II . 53 

20  Shock  Distance  Versus  Time  in  MISERS  BLUFF  II .  54 

21  Locations  of  History  Observations  in  MISERS  BLUFF  II  .  55 

22  Fitting  of  Pressure  at  89.7  m  in  MISERS  BLUFF  II . 55 

23  Fitting  of  Pressure  at  103  m  in  MISERS  BLUFF  II . 37 

24  Fitting  of  Pressure  at  114  m  in  MISERS  BLUFF  II . 5g 

25  Fitting  of  Pressure  at  120.3  m  in  MISERS  BLUFF  II . 59 


5 


LIST  OF  FIGURES  (CONTINUED) 

Figure  Page 

26  Fitting  of  Pressure  at  150.3  m  in  MISERS  BLUFF  II . 60 

27  Fitting  of  Pressure  at  174  m  MISERS  BLUFF  II . 61 

28  Fitting  of  Pressure  at  206  m  in  MISERS  BLUFF  II . 62 

29  Fitting  of  Pressure  at  211  m  in  MISERS  BLUFF  II . 63 

30  Log-log  Plots  of  History  Parameters  versus  Distance  in  MISERS 

BLUFF  II . 65 

31  Pressure  History  at  120  m  in  MISERS  BLUFF  11/(89.7-174  m) .  .  .  .  66 

32  Velocity  History  at  120  m  in  MISERS  BLUFF  II/ (89 . 7-174m)  ....  67 

33  Density  History  at  120  m  in  MISERS  BLUFF  11/(89.7-174  m)  .  .  .  .  68 

34  Dynamic  Pressure  History  at  120  m  in  MISERS  BLUFF  II/ 

(89. 7-. 74  m) . 69 

35  Pressure  History  at  120  m  in  MISERS  BLUFF  11/(103-174  m)  .  .  .  .  70 

36  Velocity  History  at  120  m  in  MISERS  BLUFF  11/(103-174  m)  .  .  .  .  71 

37  Density  History  at  120  m  in  MISERS  BLUFF  11/(103-174  m) .  72 

38  Dynamic  Pressure  History  at  120  m  in  MISERS  BLUFF  II/ 

(103-174  m) . 73 

39  Pressure  History  at  150  m  in  MISERS  BLUFF  11/(120.3-206  m)  .  .  .  74 

40  Velocity  History  at  150  m  in  MISERS  BLUFF  11/(120.3-206  m)  .  .  .  75 

41  Density  History  at  150  m  in  MISERS  BLUFF  11/(120.3-206  m) .  .  .  .  75 

42  Dynamic  Pressure  History  at  150  m  in  MISERS  BLUFF  II/ 

(120.3-206  m) . 77 

43  Estimated  Positive  Overpressure  Duration  in  a  0.1  kton  TNT 

Experiment . 79 

44  Range  Intervals  for  a  0.1  kton  TNT  Experiment . 


6 


1.  INTRODUCTION 


For  experimental  studies  of  target  response  to  high  energy  blast,  one 
needs  an  accurate  definition  of  the  blast  field  which  provides  the  load  on  the 
target.  The  general  features  of  the  flow  in  a  blast  field  can  be  predicted  by 
solving  numerically  the  governing  equations  of  the  flow.  However,  the  accuracy 
of  such  calculations  is  limited  because  the  early  stages  of  an  explosion  are 
difficult  to  model,  and  the  results  of  the  early  stage  calculations  influence 
subsequent  computations  for  later  times.  Therefore,  direct  measurements  of 
the  complete  flow  field  would  be  most  desirable;  but  such  measurements  are 
usually  restricted,  for  technical  reasons,  to  pressure  history  observations, 
and  to  shock  arrival  time  and  incident  shock  pressure  measurements  at  various 
stations.  Hence,  one  has  the  task  to  compute  other  flow  variables;  e.g.,  the 
density  and  the  particle  velocity,  from  the  measured  pressures. 

Directly  computable  from  pressure  history  observations  is  the  history  of 
the  dynamic  pressure  (internal  energy  density) ,  if  simultaneous  measurements 
of  incident  and  stagnation  pressures  are  available.  However,  the  results  are 
very  sensitive  to  observational  inaccuracies.  Typically,  the  relative  standard 
error  of  the  calculated  dynamic  pressure  is  ten  times  as  large  as  the  relative 
standard  errors  of  the  pressure  observations^.  Also,  the  calculations  do  not 
provide  any  other  flow  variables  besides  the  dynamic  pressure. 

Generally,  one  can  improve  the  accuracy  of  computed  results  by  using 
additional  information.  For  the  present  problem,  such  information  is  available 
from  two  sources:  from  pressure  history  observations  at  stations  adjacent  to 
the  station  of  interest  and  from  the  governing  equations  of  the  flow.  Compared 
to  the  governing  equations  of  the  early  stages  of  the  explosion,  the  governing 
equations  for  late  stage  flow  are  relatively  simple  and  amenable  to  numerical 
solution.  Times  and  distances  that  are  of  interest  for  target  response  studies 
typically  are  within  this  flow  stage.  Therefore,  one  can  formulate  the  prob¬ 
lem  as  a  task  to  solve  numerically  the  governing  equations  with  boundary  con¬ 
ditions  derived  from  pressure  history  and  shock  observations. 

In  this  formulation,  the  task  is  a  mathematically  ill-posed  problem  because 
the  boundary  conditions  overdetermine  the  solution  in  some  parts  of  the  flow 
field,  and  at  the  same  time  may  not  be  sufficient  to  compute  the  complete  flow 
history  for  the  full  duration  of  a  pressure  history  observation  at  some  other 
station. 

Regularizations  of  the  problem  have  been  achieved  by  using  two  opposite 
approaches.  In  one  approach,  one  deletes  all  data  that  would  overdetermine 
the  problem.  In  the  other  approach,  one  uses  all  pertinent  data,  but  deletes 
one  of  the  governing  equations,  using  it  later  for  control  calculations.  The 
former  approach  has  been  used  by  Gottlieb  and  Ritzel^.  The  latter  approach  is 
described  in  this  report. 


1 

George  D .  Tee'll  "Free-Field  Airblast  Definition  -  Event  DICE  THROW,  tT  Proceedings 
of  the  DICE  THROW  Symposium  21-23  June  1977 ,  Volume  I,  Defence  Nuclear  Agency 
Report  DNA  4377P-1 ,  pp.  7-76 ,  July  1977. 
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u James  J.  Gottlieb  and  David  V .  Ritzel ,  " Flow  Properties  of  a  Spherical  Blast 
Wave, "  presentation  at  the  Sixieme  Symposium  International  sur  les  applications 
militaires  de  la  simulation  de  souffle ,  Cahors,  France,  25-27  June  1979. 
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An  advantage  of  the  approach  by  Gottlieb  and  Ritzel  is  that  the  problem 
is  immediately  reduced  to  a  standard  task  in  numerical  mathematics.  A  disad¬ 
vantage  is  that,  in  a  simplistic  application  of  the  method,  all  observational 
inaccuracies  enter  directly  into  the  calculations  and  might  be  amplified.  In 
order  to  obtain  reasonable  results,  one  has  to  smooth  judiciously  the  data  or 
to  introduce  viscosity  into  the  computing  scheme.  A  smoothing  of  data  can  be 
done  in  different  manners  so  that  in  effect  one  may  have  to  manipulate  the 
data  in  order  to  obtain  reasonable  results.  The  data  that  are  not  used  for 
smoothing  or  as  boundary  data,  could  be  used  for  a  subsequent  check  on  the 
calculations,  but  results  of  such  tests  have  not  been  reported  in  Reference  2. 

In  our  approach,  the  problem  is  regularized  by  discarding  one  of  the 
governing  equations,  but  using  all  pressure  observations  within  a  region  of 
interest.  We  assume  that  the  observations  suffice  to  determine  a  pressure 
field  function  p(r,t)  within  the  region.  The  function  is  found  by  a  least 
squares  model  fitting,  which  takes  the  role  of  the  data  smoothing  in  the  Gottlieb 
and  Ritzel  approach.  The  pressure  field  function  then  is  substituted  into  the 
governing  equations  which  in  turn  determine  the  other  flow  variables.  Problems 
of  this  type  were  considered  by  Makino^  who  observed  that  one  does  not  need 
the  continuity  equation  for  the  flow  calculation  if  p(r,t)  is  known.  We  have, 
in  essence,  followed  Makino!s  theoretical  ideas  and  established  computer  pro¬ 
grams  that  compute  the  flow  in  the  aforementioned  manner.  The  continuity  equa¬ 
tion  is  used  at  the  end  of  the  calculations  to  check  the  accuracy  of  the  results. 
We  also  carried  out  an  analysis  of  the  sensitivity  of  the  results  to  observational 
inaccuracies,  and  included  computations  of  accuracy  estimates  in  the  computer 
programs.  In  typical  examples,  the  results  were  found  to  be  quite  insensitive 
to  this  type  of  inaccuracy.  For  example,  the  relative  standard  error  of  the 
dynamic  pressure  was  found  to  be  of  the  same  order  or  less  than  the  relative 
standard  error  of  the  overpressure  observations. 

This  report  outlines  the  theory  of  the  flow  calculations  from  pressure 
measurements  (Sections  2  through  6)  and  demonstrates  its  application  in  two 
examples  (Sections  7  and  8) .  The  first  example  treats  a  theoretically  computed 
strong  blast  field.  The  second  example  is  an  application  of  the  technique  to 
real  observations. 

A  description  of  the  use  of  the  computer  programs  will  be  published  in  a 
forthcoming  users1  manual. 


2.  BASIC  ASSUMPTIONS  AND  THEORY 

We  seek  to  determine  certain  parts  of  the  flow  field  within  a  blast  bubble 
in  air.  The  area  of  interest  is  a  relatively  narrow  strip  in  the  r,t-plane 
behind  the  initial  shock  trajectory  at  a  distance  where  the  shock  strength  is 
only  moderate.  We  shall  assume  that  the  following  conditions  are  satisfied 
within  the  area  of  interest: 

(A)  the  flowing  medium  is  an  ideal  gas  with  zero  viscosity  and  no  heat 
conduction,  and 


7 

Ray  C.  MakinOj  " An  Approximation  Method  in  Blast  Calculations , "  BRL  Memorandum 
Report  1034 3  February  1956.  (AD  #114875) 
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(B)  the  event  is  spherically  symmetric  and  the  flow  has  only  a  radial 
velocity  component  u. 

The  first  assumption  is  satisfied  in  most  applications  because  typically 
the  maximum  overpressure  at  the  target  is  only  of  the  order  of  one  megapascal. 
Within  this  pressure  regime,  air  behaves  like  an  ideal  gas.  The  second  condi¬ 
tion  is  nearly  satisfied  in  most  experiments,  because  usually  the  explosion 
source  and  the  targets  are  positioned  on  the  same  plane,  and  the  blast  bubble 
is  a  hemisphere.  Deviations  from  spherical  flow  symmetry  within  the  bubble 
may  be  caused  by  local  surface  disturbances,  by  wind,  and  by  the  presence  of 
dust  in  the  flow  near  the  ground  surface.  The  present  technique  cannot  be 
applied  to  cases  where  such  disturbances  are  not  negligible. 

The  governing  equations  for  a  flow  satisfying  the  conditions  (A)  areS 

the  continuity  equation 


—  +  p  div  u  =  0  , 
the  momentum  equation 
P  +  grad  p  =  0  , 


the  energy  equation 


P 


de  p  dp  _  » 
dt  "  p  dt 
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and  the  equation  of  state 

e  -  . 

Y-l  P 


(2.1) 


(2.2) 


(2.3) 


(2.4) 


In  Eqs.  (2.1)  through  (2.3)  we  have  used  the  material  derivative  symbol 


+  (u-grad)  .  (2.5) 

Because  the  flow  is  assumed  to  be  spherically  symmetric  (condition  (B)),  we 
also  have  the  definitions 


^ Richard  von  Mises,  " Mathematical  Theory  of  Compressible  Fluid  Flow,  ’’  Academic 
Press 3  NY,  1958. 
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(2.6) 


div  u 


1_ 

2 

r 


and 


grad  p  = 


9r  ' 


(2.7) 


Eliminating  the  specific  internal  energy  e  between  Eqs.  (2.3)  and  (2.4) 
one  obtains 


1  dp  1  dp 

pdt'YpdT 


(2.8) 


Eq.  (2.8)  can  be  integrated  along  a  particle  path  line.  The  result  is  the  well 
known  formula  for  a  particle  in  an  adiabatic  flow: 


(2.9) 


where  the  subscript  A  indicates  reference  values  at  a  point  A  on  the  particle 
path. 


The  momentum  Eq.  (2.2)  can  be  reformulated  by  substituting  in  it  the 
expression  (2.9).  The  result  is 


du 

dt 


?r 


9£ 

9r 


(2.10) 


If  the  pressure  function  p(r,t)  is  given,  e.g.,  by  measurements,  then  the 
right-hand  side  of  Eq.  (2.10)  is  a  known  function  of  r  and  t.  In  that  case, 
the  equation  can  be  integrated  numerically  together  with  the  path  line  equation 


dr  = 
dt 


(2.11) 


The  integration  provides  the  path  line  starting  at  a  point  A  and  the  particle 
velocity  along  it.  The  density  along  the  same  path  line  is  given  by  Eq.  (2.9). 
All  other  flow  variables,  such  as,  internal  energy,  dynamic  pressure,  and 
sound  speed  can  be  computed  from  p,  u,  and  p. 
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The  continuity  Eq.  (2.1)  is  not  needed  for  the  described  calculation  of 
the  flow  corresponding  to  an  observed  pressure  field  p(r,t).  One  can  use  the 
equation  to  test  the  calculated  results,  as  suggested  by  Makino^.  In  fact,  if 
the  pressure  p(r,t)  is  measured  precisely  then  this  test  provides  a  check  of 
the  validity  of  the  assumptions  (A)  and  (B)  about  the  flow  field.  In  practice, 
test  calculations  based  on  the  continuity  equation  cannot  provide  exactly  the 
same  result  as  the  integration  along  path  lines  because  the  pressure  field 
function  p(r,t)  contains  observational  as  well  as  systematic  errors.  The 
effects  of  the  former  are  estimated  in  our  approach  using  input  information 
about  the  data  accuracy.  Systematic  errors  may  manifest  themselves  by  differ¬ 
ences  between  original  and  control  calculations  that  are  larger  than  predicted 
by  the  observational  error  analysis. 

A  control  calculation  based  on  the  continuity  equation  can  be  carried  out 
as  follows.  First,  we  use  Eq.  (2.8)  and  reformulate  the  continuity  Eq.  (2.1), 
obtaining 


div  u  +  —  ”^7  =  0 
YP  dt 


(2.12) 


or 


9  /  2  v  ,  2  v 

(r  u)  4-  (r  u) 


i_  12.  +  ll  =  o 

yp  9r  YP  3t 


(2.13) 


Eq.  (2.13)  expresses  the  dependence  of  the 
A  formal  integration  of  the  equation  along 


quantity  r  u  on  r  for  t  -  const, 
a  line  t  «  const,  yields 


vl/Y 


2  . 
r  u(r,t)  = 


P(r,t) 


>  ,  i  r 

:uc  +  y  J 


(  \1/Y 

.2  I  p(g,t)  \  1  3p(£,t) 

’  \  Vn  /  p(5.t)  3t 


d£ 


(2.14) 


or 


u(r,t)  =  u  I  — 


■A?)‘ (ifcf  * 


r  ■ 

- 1 _ ^  /Y-p« 

2  ,  ,.i/y  /  p  ’  at 

r  y  p(r,t)  J 


(2.15) 


The  subscript  C  in  Eqs.  (2.14)  and  (2.15) 
C  with  the  coordinates  (rQ,t).  Using  Eq. 
velocity  u(r,t)  by  a  numerical  quadrature 
value  u  =  u(r  , t)  and  the  pressure  field 


indicates  function  values  at  a  point 
(2.15)  one  can  calculate  the  particle 
along  t  =  const. ,  if  an  initial 
function  p(r,t)  are  known. 


ZIbid. 
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In  summary,  we  proceed  as  follows  for  the  calculation  of  the  flow  field. 
First,  we  establish  a  pressure  field  function  p(r,t)  by  data  fitting,  whereby 
all  data  are  used  that  are  located  in  the  vicinity  of  the  point  of  interest 
say  B.  Next,  Eqs.  (2.10)  and  (2.11)  are  integrated  along  a  particle  path  AB, 
as  shown  in  Figure  1.  The  integration  produces  the  velocity  ug  at  B.  The 
density  pB  can  be  computed  using  Eq.  (2.9),  once  the  path  line  is  established. 
(The  flow  variables  uA  and  pA  on  the  shock  are  known  from  the  pressure  field 
function  and  shock  relations.)  Finally  the  calculated  velocity  uR  is  compared 
with  another  calculation  using  Eq.  (2.15),  applied  along  the  line  BC.  The 
velocity  u^  at  the  point  C  is  again  obtained  from  shock  relations. 


3.  NUMERICAL  INTEGRATION  AND  ACCURACY  ESTIMATES 

In  most  applications,  one  needs  the  flow  history  at  some  fixed  distance, 
saY  rg-  We  obtain  the  history,  i.e.,  the  values  of  flow  variables  at  a  series 
of  points  along  the  line  r  =  rB,  by  integrating  Eqs.  (2.10)  and  (2.11)  along 
a  number  of  path  lines,  each  starting  at  a  different  point  of  the  shock.  The 
test  calculation  of  the  velocity  is  done  by  integration  of  Eq.  (2.15)  along 
appropriate  lines  t  =  const.  Figure  2  shows  schematically  the  integration 
lines  and  the  location  of  the  computed  nodes  in  the  r,t-plane.  The  values  of 
the  flow  variables  at  the  shock  as  well  as  the  pressure  field  function  behind 
the  shock  that  are  needed  for  these  integrations,  are  obtained  by  model  fitting 
of  shock  and  pressure  observations,  respectively. 

Let  the  result  of  the  shock  model  fitting  be  certain  functions  of  the 
radial  distance  r  providing  shock  arrival  time  t  ,  shock  overpressure  p  ,  shock 
density  p^,  and  particle  velocity  u  :  S  s 


t 

s 


tg(r;e)  , 


Ps  =  Ps(r;0) 

ps  =  ps(r;G) 


(3.1) 


u  =  u  (r ; 6)  . 

s  s  ’ 


In  these  functions,  0  is  a  model  parameter  vector.  We  shall  give,  in  Section 

4,  explicit  expressions  for  the  functions  t  through  u  . 

s  s 

The  model  fitting  of  the  observed  pressure  histories  produces  a  pressure 
field  function 


p  =  pf(r,t;0)  , 


(3.2) 


Figure  1.  Integration  Paths 


The  flow  variables  are  obtained  at  the  Point  B  by  integration 
along  the  path  line  AB.  The  computed  velocity  is  checked  at 
B  by  an  integration  along  the  line  CB. 
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A  t 


Figure  2.  Computation  of  Flow  History  at  a  Given  Distance 


The  pressure  field  function  is  determined  within  the  indicated  domain 
from  pressure  history  measurements  along  the  lines  AA3,  BB3,  and  CC3, 
and  from  shock  observations.  For  r  =  rg,  the  flow  history  can  be 
calculated  between  B  and  B2,  and  test  calculations  can  be  carried  out 
between  B  and  B-j . 
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that  describes  the  overpressure  within  a  limited  region  behind  the  shock.  This 
function  will  be  handled  in  detail  in  Section  6.  The  pressure  field  function 
and  the  shock  pressure  function  are  related  by 


pg(r;9)  =  pf (r ,tg(r ,0) ;0) 


(3.3) 


The  differential  equations  for  the  path  line,  Eqs.  (2.10)  and  (2.11),  are 
in  terms  of  these  functions: 


dr 

dt 

du 

dt 


=  u 


*  F(r,t;0)  , 


(3.4) 


where 


F (r ,t;0)  = 


'PAI/V 


= 

9r 


"  k  (p ) 

_  1  X 

Ps^rA;0)  \pf »tJ0)+Po/ 


8pf (r,t;0) 
9r 


(3.5) 


and  p^  is  the  ambient  pressure. 

The  numerical  integration  of  the  system  (3.4)  is  done  by  the  following 
fourth  order  algorithm: 


tk+l 

fck 

+ 

At, 

rk+l 

rk 

+ 

At  • 

\  + 

At2  ti^/2  + 

Uk+1  = 

Uk 

+ 

At  • 

(uk  + 

-W/2  • 

Uk+1  = 

3 

at 

(W 

+  Vi 

.  h  <W 

Uk+1  = 

A 

Uk+1 

+  At 

2  ... 

(\  - 

"V,1!  - 

3 


(3.6) 


The  time  derivatives 


u,  and  u* 
k  k 


are  defined  for  this  algorithm  by 
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uk  =  F(rk’tk;0)  ’ 


(3.7) 


\ 


u,  = 


3_ 

3t 


<v 


+  U, 


The  integration  errors  of  the  algorithm  can  be  reduced  to  desired  levels 
by  monitoring  the  time  steps  At.  However,  the  accuracy  of  the  solution  is 
also  influenced  by  the  accuracy  of  the  pressure  observations  which  are  used  to 
establish  the  pressure  functions  pf(r,t;0)  and  ps(r;0).  This  sets  a  limit  on 
the  achievable  accuracy  of  any  solution  of  Eq.  (3.4). 

The  basis  for  the  accuracy  estimate  is  an  estimate  of  the  variances  (or 
standard  errors)  of  the  overpressure  observations,  which  we  assume  to  be 
obtained  by  an  analysis  of  the  measurement  process.  The  least  squares  model 
fitting  routines^  use  the  data  variances  to  calculate  an  estimate  of  the 
variance-covariance  matrix  Vq  of  the  parameter  vector  0.  Using  Vq ,  estimates 
of  standard  errors  of  pf(r,t;0)  and  pg(r;0)  can  he  obtained  by  applying  the 
linearized  law  of  variance  propagation  to  the  formulas  for  pf  and  ps,  respec¬ 
tively.  Thus,  an  estimate  of  the  standard  error  e  of  pf  is 


e 

P 


1/2 


(3.8) 


where  3pf/36  is  the  Jacobian  matrix  of  the  function  Pf(r,t;0).  (The  Jacobian 
matrix  is  a  vector  in  this  case  because  p^  is  a  scalar.) 

In  order  to  calculate  an  estimate  of  the  standard  error  of  the  velocity  u 
by  a  corresponding  formula  one  needs  to  know  the  derivative  vector  3u/30. 

Unlike  3pf/30,  that  vector  cannot  be  obtained  by  formal  differentiation  because 
u  is  not  given  by  a  formula  but  obtained  by  solving  numerically  the  equation 
system  (3.4).  Therefore,  we  differentiate  that  system  producing  another  system 
of  differential  equations  where  the  unknown  functions  are  the  derivatives 
Bu/30  and  3r/30.  The  new  system  is 


d  /3r\  =  3u  \ 

d  t  \3  0  /  30  ’  1 

[  (3.9) 

d  /3u\  3F  3F  3r  | 

dt  \30  /  30  +  3r  30  ’  J 


where 


Aivars  Celmiys,  "A  Manual  for  General  Least  Squares  Model  Fitting USA 
ARRADCOM /Ballistic  Research  Laboratory  Technical  Report  ARBRL-TR-02167 ,  June 
1979.  (AD  #B040229L) 
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3F  3F  3r 
30  3r  30 


3P, 


3P, 


+ 


30  '  Y(PA+P0)  36 


1 

YCPf+Po) 


3p  3p 

W  +  "17 


+ 


3pf/3r 


~2 

3  pf 
3r  30 


3jr 

30 


(3.10) 


Eqs.  (3.9)  are  integrated  numerically  concurrently  with  the  path  line 
Eqs.  (3.4),  and  using  the  same  time  steps  At.  For  simplicity,  and  because  high 
accuracy  is  not  needed,  we  use  a  third-order  scheme  for  the  integration  of 
Eqs.  (3.9).  The  third-order  scheme  is  obtained  from  the  algorithm  (3.6)  by 
setting  ti^  and  ti^.  equal  to  zero. 

The  initial  values  for  the  integration  of  Eqs.  (3.9)  are 


\30/r=r 

A 


0 


9 


and 


3_ 

30 


Us(rA;e) 


(3.11) 


The  derivatives  3r/30  are  zero  at  the  starting  point  of  the  path  line  because 
r=r  is  a  prescribed  initial  value. 

The  integration  of  Eq .  (3.10)  with  the  initial  conditions  (3.11)  yields 
numerical  values  of  the  vectors  3u/30  and  3r/30  at  any  point  along  the  path 
line.  Of  particular  interest  is  the  accuracy  of  the  flow  at  the  end  point, 
where  the  line  intersects  with  r=rg.  The  end  point  itself  has  an  uncertainty 
in  the  t-direction  only  that  can  be  estimated  by  the  linearized  law  of  vari¬ 
ance  propagation,  if  the  derivative  3t/30  is  known.  We  obtain  a  formula  for 
3t/30  at  the  end  point  as  follows.  For  a  given  path  line  the  time  value  t^  at 

r  is  a  function  of  the  type 
B 


CB  ■  CA  +  T(rA-rB> 


(3.12) 


Its  derivative  with  respect  to  0  is 

9tB  =  _^A  dT_  /  3r\  3T_  /3r^ \ 

30  30  +  3r  \ 30 /r=r  Sr,.  \30/r=r„ 

A  N  A  B  B 


or,  because  of  Eq.  (3.11), 


30 


_ 3_ 

30 


B 


(3.13) 
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The  derivatives  9pg/90,  that  are  needed  for  the  estimation  of  the 
accuracy  of  p  at  the  end  point,  can  be  calculated  by  a  differentiation  of  Eq. 
(2.9):  " 


!!b  _ 

90  PB 


■fe 


90 


9p 


A+  1 


8Pr 


Y(pA+P0)  86  y(pB+Po^  80 


(3.14) 


where 


90  90  s  v  A 


ps(^a;0)  . 


(3.15) 


9p 


A  9  , 

=  To-  P.(rAJ0) 


90  90  *s  A’ 


(3.15) 


and 


3Pj 

90" 


90  Pf^rB,t:B;0^  +  9t  [Pf  ^rB’  tB’9^] 


90 


(3.16) 


The  flow  at  the  end  point  of  the  path  line,  i.e.,  at  its  intersection 
with  the  line  r=r  is  defined  by  the  four  component  vector 

D 


H  =  (V  PB ’  V 


(3.17) 


The  accuracy  of  the  vector  can  be  characterized  by  a  variance-covariance  matrix 

V  of  the  components  of  H,  that  can  be  calculated  from  V  by  the  formula 
H  0 


Vfl  estimates  the  uncertainties  of  the  flow  field  that  are  caused  by  measurement 
inaccuracies.  The  elements  of  the  Jacobian  matrix  9H/9<f)  in  (3.18)  are  given 
by  the  solution  of  the  equation  system  (3.11)  and  the  formulas  (3.13),  (3.14), 
and  (3.16),  respectively. 

4.  DETERMINATION  OF  SHOCK  FUNCTIONS  BY  MODEL  FITTING 

In  order  to  compute  the  complete  blast  field  by  the  method  described  in 
previous  sections,  one  needs  the  pressure  field  function  pf(r,t;0)  within  the 
region  of  interest.  We  determine  that  function  by  a  model  fitting  to  pressure 
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observations,  most  of  which  consist  of  overpressure  history  recordings  at  di¬ 
verse  stations.  The  initial  value  of  most  overpressure  history  recordings  is 
the  incident  shock  overpressure.  From  the  viewpoint  of  data  handling,  it  is 
advantageous  to  use  these  initial  values  separately  from  the  rest  of  pressure 
data  for  a  determination  of  the  shock  functions  (3.1).  Then  one  can  obtain  a 
shock  description  that  is  accurate  within  a  large  distance  range  and  use  the 
same  shock  functions  for  different  field  fittings,  as  will  be  described  in  the 
next  sections.  Also,  in  a  separate  shock  fitting  one  can  easily  make  use  of 
incomplete  records  from  which  either  the  pressure  observation  or  the  shock 
arrival  time  observation  is  missing.  (Incomplete  records  do  occur  occasionally 
in  blast  field  experiments.)  The  use  of  all  available  data  makes  the  shock 
functions  less  sensitive  to  individual  observational  inaccuracies.  Since  the 
shock  functions  also  can  be  approximated  very  accurately  by  simple  expressions, 
a  separate  shock  fitting  is  more  practical  than  a  concurrent  determination  of 
shock  and  field  functions.  Details  of  the  shock  fitting  algorithm  will  be 
given  in  this  section.  The  complementary  pressure  field  fitting  will  be  de¬ 
scribed  in  Sections  5  and  6. 

The  medium  into  which  the  shock  propagates  is  assumed  to  be  an  ideal  gas, 
characterized  by  its  molar  mass  M  and  the  ratio  y  of  its  specific  heats.  We 
describe  the  initial  state  of  the  gas  by  its  pressure  p  and  temperature  T  . 

The  corresponding  initial  density  is 


(4.1) 


P  o  T  *  R 

°  g 


where  Rg  =  8.31434  J/(mole*K)  is  the  universal  gas  constant.  The  sound  speed 
in  the  ambient  gas  is 


(4.2) 


Let  pg  be  the  incident  shock  overpressure.  Then  the  shock  velocity, 
particle  velocity  and  density  behind  the  shock  are  given  in  terms  of  p  by  the 
following  set  of  formulas^. 


Shock  velocity: 


U  =  c 


o 


I±l!£/1+^x!£\l1/2 


2  po  V  v+1  Ps ) 


(4.3) 


Particle  velocity  behind  the  stock: 


2 


(4.4) 


4 


Ibid.  p.  200  ff. 
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Density  behind  the  shock: 


P 


s 


Y+l  Ps  \ 

2Y  pJ 


iziZnY1 

P„/ 


y+l 

p  - - 

o  y-1 


2y  Po\ 

Y+1  ps  / 


2y_  £o  \  1 
Y-1  Ps  )  ' 


(4.5) 


The  shock  arrival  time  tg(r)  at  a  distance  r  from  the  center  of  the  ex¬ 
plosion  and  the  shock  velocity  U  are  related  by 


t 

s 


(r) 


(4.6) 


For  the  numerical  treatment  of  shock  data  we  choose  a  rational  function 
of  r  as  an  approximation  of  the  shock  overpressure  function  pg(r).  The  rational 
function  has  three  parameters  and  an  asymptotic  behavior  for  small  and  large 
r  as  predicted  by  blast  theory.  The  function  is 


Pg(r ;a,b,c) 


=  a/r  +  b/r2  +  c/r3 


(4.7) 


The  corresponding  shock  arrival  time  function  ts(r)  is 
tuting  the  expression  (4.7)  into  the  shock  velocity  formula 
suiting  function  into  the  integral  (4.6).  The  final  result 


obtained  by  substi- 
(4.3)  and  the  re- 
is 


tg(r ;a,b,c,d) 


dx 


c 

o 


Vl  +  £^-(a/x  +  b/x2 

T  2YPo 


+  c/x3) 


(4.8) 


The  integral  in  Eq.  (4.8)  must  be  evaluated  numerically,  and  the  evaluation 
provides  the  shock  arrival  time  t  at  any  given  distance  r.  (In  the  computer 
programs  for  shock  fitting,  the  quadrature  is  done  by  a  Romberg  algorithm.) 

The  lower  limit  r^  of  the  integral  in  Eq.  (4.8)  can  be  chosen  arbitrarily, 
e.g.,  it  may  be  set  equal  to  the  smallest  observed  distance.  Once  r^  is  fixed, 
then  both  model  functions,  Ps(r)  and  tg(r),  are  completely  determined  by  the 
four  components  of  the  model  parameter  vector  0  =  (a,b,c,d)T. 

Appropriate  values  of  the  four  parameters  were  found  by  a  least  squares 
data  fitting  that  uses  the  shock  observations  as  data.  Each  shock  observation 
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is  a  three  component  vector  (Ps,  Ts,  Rs)  ,  i.e.,  the  observed  shock  overpressure, 
arrival  time,  and  distance.  The  model  constraints  for  the  observations  are 
Eqs.  (4.7)  and  (4.8).  They  may  be  formulated  as  the  equation  system 


(4.9) 


or  in  some  other  equivalent  form  consisting  of  two  simultaneous  equations  for 


each  three-component  vector  of  observations.  Problems  of  this  type  can  be 
solved  numerically  with  the  aid  of  the  utility  program  COLSMU^,  that  was  also 
used  in  the  present  task.  The  program  is  sufficiently  flexible  to  accommodate 
observation  sets  from  which  either  the  pressure  or  the  time  observation  is 
missing.  (The  constraint  for  such  an  incomplete  set  consists  of  only  one  of 
the  two  equations  (4.9).) 

The  first  constraint  in  Eq.  (4.9)  was  formulated  for  the  model  fitting 
task  as  follows: 


3  2 

f1(p,t,r;  a,b,c,d)  =rp-ra-rb-c 


0  . 


(4.10) 


The  first  order  derivatives  of  f^  are 


3f1/9(p,t,r)  =  (r3,  0,  3r2p  -  2ar  -  b) 


(4.11) 


and 


8f1/8(a,b,c,d)  =  (-r2,  -r,  -1,  0)  . 


(4.12) 


The  non-zero  second  order  derivatives  of  f^  are 


92f1/(8r8p)  =  3r2  , 


(4.13) 


92f1/(9r9b)  =  -1  . 


5 


Ibid . 
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The  second  constraint  in  Eq.  (4.9)  was  used  in  the  model  fitting  task  in 
the  following  form 


f9(p»t,r;  a,b,c,d)  = 


7 


dx 


v 


1  +  T* (a/x+b/x  +c/x  ) 


- +  (d-t)c  =  0  ,  (4.14) 

2,  :  3,  ° 


where 


T  =  (y  +  l)/(2ypo)  .  (4.15) 

Let 

Q(x)  =  ^ 1  +  T  • (a/x  +  b/x2  +  c/x3)  .  (4.16) 

Then  the  first  order  derivatives  of  the  model  function  f^  are 

9f2/3(p,t,r)  =  (0,  -cq,  1/Q(r))  (4.17) 


The  non-zero  second  order  derivatives  of 


f2  are 


92f2/3r2  =  j  T  *(a+2b/r+3c/r2)/(r2Q3(r))  , 


(4.19) 


32f2/(3r9a)  =  -  j  T/(rQ3(r))  , 

92f2/(9r3b)  =  -  ^  r/(r2Q3(r))  ,  l  (4.20) 

32f2/(3r3c)  =  -  |  T/(r3Q3(r))  , 
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(4.21) 


2 


/ (3b3c)  =  7-  T2  f< 


dx/ (x5Q5(x)) 


92f2/9c2  =  f  r2 


The  limits  of  the  integrals  in  Eq.  (4.21)  are  the  same  as  in  the  constraint 
Eq.  (4.14).  All  integrals  are  evaluated  numerically  by  Romberg  algorithms. 


5.  ANALYSIS  OF  PRESSURE  HISTORIES 


The  data  for  a  determination  of  the  pressure  field  function  consist  of  a 
number  of  pressure  histories  observed  at  various  locations  within  an  area  of 
interest.  Before  determining  a  field  function  pf(r,t)  that  describes  the  pres¬ 
sure  p  as  a  function  of  r  and  t  within  the  whole  area,  we  find  functions  p^(t) 
that  fit  the  observed  individual  pressure  histories.  The  purpose  of  the  fitting 
of  the  individual  histories  is  to  obtain  trends  of  the  histories  for  increasing 
distance  from  the  explosion.  These  trends  are  needed  to  construct  the  pressure 
field  function  pf(r,t),  as  will  be  described  in  the  next  section. 

The  model  function  chosen  to  fit  the  individual  overpressure  histories  is 


(5.1) 


where  A,  B,  and  C  are  model  parameters,  pg  is  the  shock  overpressure,  and  tg 
is  the  shock  arrival  time.  The  values  of  pg  and  tg  depend  on  the  location  of 
the  pressure  gage,  i.e.,  on  the  distance  r.  They  are  calculated  using  the 
formulas  (4.7)  and  (4.8). 

The  data  consist  of  two-component  vectors  (p,t)  giving  the  overpressure  p 
for  various  time  values  t.  The  model  fitting  problem  can  be  solved  numerically 
by  the  utility  routine  COLSAC5.  The  resulting  parameters  A,  B,  and  C  are  valid 
only  for  the  particular  history  used  as  input,  and  they  have  no  general  physi¬ 
cal  significance. 


5 


Ibid. 


23 


The  constraint  function  for  the  model  fitting  of  history  data  was  formu¬ 
lated  as  follows 

f(p,t;  A,B,C)  =  (Pg-C) » exp  (a*  (t-tg)  +  B«(t-tg)2)  +  C  -  p  .  (5.2) 

Let 

E  =  exp  (A*(t-tg)  +  B«(t-tg)2)  .  (5.3) 

Then  the  first  order  derivatives  of  the  constraint  function  (5.2)  are 

9f/9(p,t)  =  ^-1,  [a  +  2B(tg-t))*  (pg-C)  •  (5.4) 

and 

9f/9(A,B,C)  =  f(t-t  )(p  -C)E,  (t-t  )2(p  -C)E,  -E  +  l)  .  (5.5) 

\  s  S  S  S  J 

The  non-zero  second  order  derivatives  are 

92f/9t2  =  [(A  +  2B(t-t  ))2  +  2b] (p  -C)E  ,  (5.6) 

s 

92f/9A2  =  (t-t  )2(p  -C)E  , 
s  s 

9 2f / (9A9B)  =  (t-t  )3(p  -C)E  , 
s  s 

92f/(9A9C)  =  -  (t-tg) E  , 

92f/9B2  =  (t-t  )4(p  -C)E  , 
s  s 

9 2f / (9B9C)  =  -  (t-tg)2E  , 
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92f/(9t9A)  =  [l  +  (A  +  2B (t-t  ))(t-t  )](p  -C)E  , 
\  S  /  s  s 


92f/(9t9B)  =  [2  +  (A  +  2B (t-t  ))(t-t  )](t-t  )(p  -C)E  ,  }  (5.8) 

\  S  /  S  S  S  I 

92f/(9t9C)  =  -  (A  +  2B(t-tg))E  . 


After  the  values  of  the  parameters  A,  B,  and  C  are  obtained  for  all  his¬ 
tories  within  the  area  of  interest,  their  trends  as  functions  of  the  distance 
r  are  analyzed.  The  analysis  consists  of  finding  for  each  parameter  a  power 
function  of  r  as  an  approximation.  Let,  for  example,  Aj[  be  the  value  of  the 
parameter  A  obtained  from  a  history  fitting  at  the  distance  r^.  Then  one  can 
determine  a  function  ArnA  that  approximates  the  data  sets  ( |  Aj_ |  ,  r-jJ  .  The 
function  is  determined  by  minimizing  an  object  function  W,  that  is  defined  as 
follows 


w  =  £  (lnlAil  -  lnl 


A I  -  n  In  r.) 2 A. 2 
1  A11 


(5.9) 


The  solution  of  the  minimization  problem  is  given  by  the  formulas 

2 


D  = 


XX2l}Ailn  "[]CAi2ln  'i 


(5.10) 


A  = 


)^  \A^lnr^)  2-^  \A2ln  |  A^  |  lnrp  E(Ailnri: 


1 


'D  ,  (5.11) 


nA 


EAE<Ailn  1  Ai 1  ^  E(Ailnri)E 


(A1ln|Ai|) 


1  1 


D 


(5.12) 


The  sign  of  A  is  set  equal  to  the  sign  of  that  parameter  A.  which  corresponds 

to  the  smallest  distance  r..  1 

1 

Trends  of  the  parameters  B  and  C  are  calculated  by  identical  formulas  that 
yield  the  values  of  B,  n^,  C,  and  nQ,  respectively.  The  exponents  n^»  ng, 
and  nc,  are  used  for  the  construction  of  the  pressure  field  function  pj(r,t), 
as  described  in  the  next  section.  The  values  of  A,  B,  and  C  are  used  as  initial 
approximations  of  certain  parameters  of  that  function. 
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Numerical  experiments  with  different  data  sets  have  shown  that  the  expon¬ 
ents  n^,  ng,  and  uq  need  not  be  determined  very  accurately.  Equally  good 
approximations  of  pressure  data  can  be  obtained  using  field  functions  with 
exponents  that  vary  within  relatively  large  ranges.  The  exponents  have  no 
physical  significance,  because  the  parameters  A^,  B^,  and  depend,  not  only 
on  the  distance  r^,  but  also  on  various  other  factors  such  as  the  duration  of 
the  observed  history,  its  noise  level,  the  number  and  distribution  of  the  ob¬ 
servations,  etc. 


6.  DETERMINATION  OF  PRESSURE  FIELD  FUNCTIONS  BY  MODEL  FITTING 

A  pressure  field  function  describes  the  overpressure  in  the  blast  field 
as  a  function  of  the  distance  r  and  the  time  t.  The  function  is  obtained  by 
data  fitting  of  pressure  observations  and  it  is,  therefore,  a  valid  approximation 
to  the  actual  pressure  field  only  within  the  limited  region  from  which  the  data 
are  selected.  In  an  analysis  of  a  given  event,  one  may  use  several  pressure 
field  functions,  each  representing  the  pressure  in  a  different  region. 

In  the  present  problem,  we  have  used  the  following  five  parameter  model 
function  as  a  representation  of  the  overpressure  field: 


pf  (r,t;A1,A2,B1,B2,C1)  =  [pg(r)-C  (r  jC^)]  exp  (Q(r .tjA^A^B^B^)  4€(r;C1)  , 

(6.1) 


where 


Q  =  [t-ts(r)](A1+A2r).rnA+J-t_t-s(r)]2(Bi+B2r) 


r*B 


(6.2) 


and 


C 


=  C  -r 


nC 


(6.3) 


The  functions  p  (r)  and  t  (r)  are  the  shock  overpressure  and  shock  arrival 
time  functions,  resplctively .  Both  are  determined  by  shock  fitting  (Section  4). 
The  exponents  n^,  n^,  and  n q  are  assumed  to  be  known  from  an  analysis  of  the 
trends  of  individual  history  fittings,  as  described  in  Section  5.  The  remaining 
five  parameters,  A]_,  A2,  B]_,  B2,  and  are  determined  by  model  fitting  to  over¬ 
pressure  observations  within  regions  of  interest. 

The  data  for  the  determination  of  the  model  parameters  are  the  observed 
overpressure,  time,  and  distance,  i.e.,  the  data  are  three  component  vectors 
(p,t,r).  The  constraint  function  is  used  for  the  data  fitting  in  the  following 
formulation: 
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(6. A) 


f (p,t,r;A1,A2,B1,B2,C1)  =  (ps-C)e^  +  C-p  , 

where  Q  and  C  are  defined  by  Eqs.  (6.2)  and  (6.3),  respectively.  The  constraint 
equations  are  scalar  and,  therefore,  the  model  fitting  problem  can  be  solved 
numerically  by  the  utility  routine  COLSAC^  for  standard  least  squares  problems. 

The  first  and  second  order  derivatives  of  the  constraint  function  f  can 
be  conveniently  expressed  by  the  following  formalism.  Let  x  and  y  be  any  two 
of  the  eight  arguments  of  f,  and  let  derivatives  with  respect  to  x  and  y  be 
denoted  by  corresponding  subscripts.  Then 


£x  -  eQ[(VC>  +  <VC)  J  -  <P-C)X  ■  (6'5> 

and 

f  =  eQ[(p  -C)  (Q  +Q  Q  )+(p  -C)  Q  +(p  -C)  Q  +(p  -C)  l-(p-C)  .  (6.6) 

xy  L  rs  Xxy  xxxy  rs  xxy  rs  yxx  rs  xyJ  r  xy 


The  derivatives  of  Q,  p^(r) ,  p  and  C(r;C^)  can  be  easily  computed. 


The  pressure  field  fitting  provides  an  optimal  set  of  the  field  parameters 
A-j_ ,  A2,  B-p  B2,  and  C-^,  and  their  estimated  variances  and  covariances.  The 
overpressure  p^,  computed  with  the  formula  (6.1),  depends,  however,  also  on 
the  four  shock  parameters  a,  b,  c,  and  d  through  the  shock  functions  ps(r)  and 
tg(r).  Therefore,  when  one  estimates  the  accuracy  of  p^,  one  has  to  take  into 
account  the  variances  and  covariances  of  all  nine  parameters,  namely,  the  four 
shock  parameters  and  the  five  field  parameters.  Let  0  be  the  vector  of  all 
nine  parameters  and  let  Vq  be  the  corresponding  variance-covariance  matrix  of 
the  components  of  0.  Then  an  estimate  of  the  variance  of  p  is  the  square  of 
Eq.  (3.8),  viz. , 


v 

P 


(6.7) 


The  matrix  Vq  was  composed  as  follows  from  the  estimated  variance-covariance 
matrix  ° ^  t*ie  s^oc^  parameters,  and  from  the  corresponding  matrix 

of  the  field  parameters: 


V 


0 


shock 


0 


Vf ield 


(6.8) 


In  the  matrix  (6.8), 
have  been  neglected. 


correlations  between  shock  parameters  and  field  parameters 
In  reality  correlations  between  both  parameter  groups 


5 Ibid. 
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exist,  because  the  shock  fitting  results  were  used  to  obtain  the  field  param¬ 
eters.  The  neglect  of  the  correlations  most  likely  results  in  an  overestimate 
of  the  pressure  variance. 

The  reason  for  the  neglect  of  the  correlations  is  programming  expediency. 
The  correlation  could  be  obtained,  e.g.,  by  a  joint  adjustment  of  all  data 
(shock  data  and  pressure  history  data)  using  the  utility  program  COLSMU5.  How¬ 
ever,  this  would  require  an  unacceptable  investment  in  programming  because  of 
the  heterogeneity  of  the  data  and  constraints.  Also,  algorithmic  difficulties 
can  be  expected  for  this  problem,  because  of  the  nonlinearities  involved,  re¬ 
quiring  additional  effort  for  their . resolution. 


7.  A  THEORETICAL  TEST  EXAMPLE 

The  blast  field  computation  method  was  tested  on  two  examples:  on  a 
theoretical  strong  shock  flow  field  and  on  real  measurements.  The  former  test 
is  described  in  this  section. 

The  theoretical  flow  field  was  obtained  by  using  the  computer  programs 
described  in  Reference  6.  The  programs  are  based  on  self-similar  solutions  of 
flow  equations7 that  approximate  the  flow  of  strong  blasts.  The  theory  of 
the  approximation  is  limited  to  blasts  in  ideal  gas  that  are  generated  by  point 
explosions,  and  for  which  the  ambient  pressure  is  negligible  compared  to  the 
incident  shock  pressure.  The  latter  assumption  makes  the  numerical  results  of 
the  theory  ambiguous,  because  the  theoretically  computed  pressure  can  be  inter¬ 
preted  either  as  overpressure  or  as  pressure.  In  our  test  example,  we  chose 
the  latter  interpretation.  Consequently,  the  overpressure  ’’observations"  for 
the  example  were  obtained  by  subtracting  the  ambient  pressure  from  the  theoreti¬ 
cal  pressure  values  which  were  supplied  by  the  strong  blast  computer  programs. 

The  explosion  in  the  test  example  was  assumed  to  be  generated  by  the 
release  of  209.2-1012  J  of  energy  (equivalent  to  a  50  kton  TNT  explosion). 

Table  1  gives  the  ambient  ideal  gas  specification  and  the  assumed  locations 
and  durations  of  overpressure  observations.  The  example  was  chosen  such  that 
the  strong  blast  assumption  was  satisfied  reasonably  well.  (The  pressure  ratio 
ranges  from  190  at  the  beginning  of  the  120  m  history  reading  to  60  at  the 
beginning  of  the  170  m  history.)  The  locations  of  the  history  readings  in  the 
r,t-plane  are  shown  in  Figure  3  and  the  assumed  data  accuracies  are  listed  in 
Table  2. 


^ Ibid. 

g 

Aivars  CelmiyS,  "Strong  Blast  Wave  Computer  Programs,  "  USA  ARRADCOM /Ballistic 
Research  Laboratory  Technical  Report  ARBRL-TR-02264,  September  1980.  (AD  #A092346) 

7 

L.I.  Sedov,  "Similarity  and  Dimensional  Methods  in  Mechanics,"  Academic  Press, 

NY,  1959. 

g 

0.  Laporte  and  T.S.  Chang,  "Curved  Characteristics  Behind  Strong  Blast  Waves," 
Physics  of  Fluids,  15,  pp.  502-504,  1972. 
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TABLE  1.  PARAMETERS  OF  THE  THEORETICAL  TEST  EXAMPLE 


Explosion 

E  =  209.2  TJ 


Ambient  Gas 

Molar  mass  M  =  28.96  g/mole 

Specific  heat  ratio  y  =  1.40 

Temperature  T  =  293  K 

Pressure  p  =  101.325  kPa 
o 

Range  of  Shock  Observations 

90-170  m 


History  Observations 

Distances 
120  m 
130  m 
140  m 
150  m 
160  m 
170  m 


Times 

11.0  -  15.0  ms 
13.4  -  17.9  ms 

16.2  -  21.5  ms 

19.2  -  27.2  ms 
22.6  -  29.6  ms 

26.3  -  32.0  ms 


TABLE  2.  ASSUMED  ACCURACY  OF  THE  THEORETICAL  DATA 

Time:  +  0.1  ms 

Pressure:  +  5% 

Distance:  +  0.1  m 
Elevation:  +  0.01  m 
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Figure  3.  Location  of  History  Observations  in 
Theoretical  Test  Example 


calculated  particle  trajectories 


theoretical  particle  trajectories 


The  leftmost  particle  trajectory  limits  the  region  where 
the  complete  flow  field  can  be  computed. 


Distances  are  expressed  in  meters  and  times  are  expressed  in  seconds. 
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Shock  arrival  times  and  pressures  were  assumed  to  be  observed  between 
90  m  and  170  m  at  2.67  m  intervals.  In  order  to  test  the  program’s  capability 
to  handle  incomplete  data  sets,  it  was  assumed  that  from  some  sets  either  the 
pressure  or  the  time  observation  was  missing.  Figures  4,  5,  and  6  show  the 
shock  data  and  the  fitting  curves  with  their  three  standard  error  confidence 
limits.  The  least  squares  values  of  the  shock  fitting  parameters  are  given  in 
Table  3. 

The  figures  show  that  the  confidence  intervals  for  the  fitted  curves  become 
very  large  outside  the  range  of  observations.  The  rapid  increase  is  caused  by 
the  redundancy  of  the  parameters  a  and  b,  which  is  manifested  also  by  their 
large  standard  errors  and  the  close  correlations  between  the  three  parameters 
a,  b,  and  c.  (See  Table  3.)  The  reason  for  the  redundancy  is  the  fact  that 
the  theoretical  shock  pressure  formula  for  the  strong  shock  is  p  =  c/r^.  One 
can  expect  that  the  increase  of  the  confidence  intervals  will  be  less  pronounced 
in  cases  where  the  observations  are  not  restricted  to  either  the  strong  shock 
regime  or  to  the  weak  shock  regime  (p«a/r). 

The  next  step  towards  the  flow  field  determination  is  a  fitting  of  indi¬ 
vidual  pressure  histories,  as  described  in  Section  5.  Figures  7  and  8  show  two 
samples  of  the  results  of  the  fitting. 

Each  history  fitting  provides  values  of  three  parameters,  A,  B,  and  C, 
which  determine  the  model  function  (5.1).  Figure  9  shows  log-log  plots  of  the 
three  parameters  versus  the  distances  at  which  the  history  recordings  were 
made.  The  field  fitting  program  calculated  trends  of  the  parameters  by  fitting 
a  power  of  the  distance  to  the  curves  shown  in  Figure  9.  The  exponents  were 
calculated  using  Eq.  (5.12)  and  the  following  results  were  obtained: 

n .  =  -2.3  , 

A 

nB  =  "4'5  » 


The  corresponding  five  parameter  overpressure  field  function  is 

Pf  =  (pg(r)  -  C1/r2'9)  +  C1/r2’9  (7.1) 


with 


Q  =  (t-ts(r))(A1+A2r)/r2,3  +  (t-tg(r))2  (B  +B  r) /r4'5  ,  (7.2) 

{  ! 

where  Ps(r)  and  tg(r)  are  the  shock  overpressure  and  arrival  time,  respectively. 
The  numerical  results  of  the  pressure  field  fitting  are  listed  in  Table  4.  The 
large  parameter  error  estimates  and  close  correlations  between  parameters  again 
indicate  that  probably  some  of  the  parameters  are  redundant  for  this  example. 
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PRESSURE 


DTSTRNCE 


BLAST  FIELD  BY  SBL-ROUTINES 

CONFIDENCE  LIMITS  FOR  3,0  STANDARD  ERRORS 
WITHOUT  THE  FACTOR  ER2  =  ,017 

ADJUSTED  ARE  OBSERVATIONS  OF 
PRESSURE*  DISTANCE  AND  TTME 


Figure  4.  Shock  Overpressure  vs.  Distance  in  Theoretical  Test  Example 


6  Complete  observation  set 
X  Set  with  missing  time  observation 

f 

Distance  is  expressed  in  meters  and  pressure  is  expressed  in  pascals. 
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PRESSURE 


BLRST  FIELD  BY  SBL-R0UT INES 

CGNEIDENCE  LIMITS  FOR  3.0  STANDARD  ERRORS 
WITHOUT  THE  E ACTOR  ER2  =  .00 

ADJUSTED  ARE  OBSERVATIONS  OF 
PRESSURE. DISTANCE  ANO  TIME 


Figure  5.  Shock  Overpressure  vs.  Time  in  Theoretical  Test  Example 
Time  is  expressed  in  seconds  and  pressure  is  expressed  in  pascals. 
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TIME 


BLAST  FIELD  BY  SBL-ROUTINES 

CONFIDENCE  LI  KITS  FOR  3.0  STANDARD  ERRORS 
WITHOUT  THE  FACTOR  ER2  =  .017 

ADJUSTED  ARE  OBSERVATIONS  OF 
PRESSURE*  DISTANCE  AND  TIME 


Figure  6.  Shock  Distance  vs.  Time  in  Theoretical  Test  Example 


■  Complete  observation  set 
X  Set  with  missing  pressure  observation 


Time  is  expressed  in  seconds  and  distance  is  expressed  in  meters. 
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TABLE  3.  SHOCK  PARAMETERS  OF  THEORETICAL  EXAMPLE 


The  model  is  given  by  Eqs.  (4.7)  and  (4.8).  The  range  of 
fitting  is  90  m  through  170  m. 

Parameters  and  Standard  Errors 

a  =  (-26.3  +  659.6) *10^  Pa*m 
b  =  (  1.13  +  171.5) -109  Pa-m2 

c  =  (32.81  +  10.94) -1012  Pa*m3 
d  =  (5.380  +  0.044)  ms  =  arrival  at  90.0  m 


The  standard  error  of  weight  one,  mQ  =  0.01674,  is  not 
included  in  the  standard  errors  of  the  parameters. 


Correlation  Matrix  of  the  Parameters 


/  1.0 

-0.99682773 

0.98588603 

0.04936359  \ 

I  -0.99682773 

1.0 

-0.99598531 

-0.08560247  \ 

\  0.98588603 

-0.99598531 

1.0 

0.13819642  I 

\  0.04936359 

-0.08560247 

0.13819642 

1.0  / 

NOTE:  The  standard  error  of  weight  one,  mQ,  depends  on  the  scatter  of  data  as 

well  as  on  systematic  differences  between  the  model  and  measurement.  For  this 
example  the  data  scatter  is  zero,  because  theoretical  flow  values  were  used 
as  data.  Therefore,  the  listed  mQ  characterizes  only  the  systematic  differ¬ 
ence  between  the  pressure  model  and  the  strong  blast  pressure.  Because  this 
difference  is  small,  the  value  of  mQ  is  also  small.  The  listed  standard  errors 
of  the  parameters  are  estimates  of  the  parameter  inaccuracies  that  would  be 
caused  by  data  scatter  (standard  errors)  as  specified  in  Table  2. 
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BLRST  FIELD  BY  SBL-ROUTINES 
SBL  R=120 

FITTED  CURVE  WITH  2 .0  STRNDRRD  ERRORS 
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Figure  7.  Fitting  of  Pressure  at  120  m  in  Theoretical  Test  Example 

(The  data  accuracy  is  indicated  by  a. two  standard  error  ellipse 
around  each  observed  point.) 


Time  is  expressed  in  seconds  and  pressure  is  expressed  in  pascals. 
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PRESSURE 


BLRST  FIELD  BY  SBL-RQUTINES 
SBL  R-170 

FITTED  CURVE  WITH  2.0  STRNDRRD  ERRORS 
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Figure  8.  Fitting  of  Pressure  at  170  m  in  Theoretical  Test  Example 

(The  data  accuracy  is  indicated  by  a  two  standard  error  ellipse 
around  each  observed  point.) 


Time  is  expressed  in  seconds  and  pressure  is  expressed  in  pascals. 
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Figure  9.  Log-log  Plots  of  History  Parameters  vs.  Distance  in  Theoretical  Test  Example 


A  test  fitting,  for  which  A2=B2=0  was  assumed,  produced  indeed  an  almost  as 
good  approximation  to  the  overpressure  field  as  the  function  defined  by  Eqs. 
(7.1)  and  (7.2).  For  subsequent  calculations,  however,  the  full  parameter  set 
of  Table  4  was  used,  because  all  five  parameters  are  likely  to  be  needed  for 
real  life  observations. 

Finally,  the  results  of  the  field  fitting  and  shock  fitting  were  used  to 
compute  flow  field  histories  at  150  m  and  160  m  from  the  explosion. 


TABLE  4.  OVERPRESSURE  FIELD  PARAMETERS  OF  THE  THEORETICAL  TEST  EXAMPLE 


The  model  is  given  by  Eq.  (7.1) 

The  total  number  of  observed  (p,t,r)-sets  is  92 


Parameters  and  Standard  Errors 


Ax  =  (-3.706  +  0.918) *107  m2'3/s 
A2  =  (2.210  +  6.111) *104  m1' 3/s 

Bx  =  (1.027  +  1.231) -1014  m4,5/s 
B2  =  (-2.384  +  8.292) -1011  m3‘5/s 
C  =  (3.840  +  0.301)-1012  m2‘9-Pa 


Standard  error  of  weight  one,  mQ  =  0.1650,  is  not 
included  in  the  standard  errors  of  the  parameters. 
(See  NOTE  in  Table  3) 


Correlation 

Matrix  of  the 

Parameters 

/  1.0 

-0.98719841 

-0.93452111 

0.92814584 

-0.09747745  \ 

/  -0.98719841 

1.0 

0.94741921 

-0.94678284 

-0.05480633  \ 

1  -0.93452111 

0.94741921 

1.0 

-0.99916227 

-0.07989523  1 

\  0.92814584 

-0.94678284 

-0.99916227 

1.0 

0.11085262  / 

'  -0.09747745 

-0.05480633 

-0.07989523 

0.11085262 

1.0 
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The  flow  histories  of  150  m  are  shown  in  Figures  10  through  14.  The  length 
ox  the  computable  history  is  limited  to  the  end  time  of  about  22  ms,  at  which 
time  the  leftmost  particle  path  reaches  the  distance  150  m  (see  Figure  3).  If 
the  field  history  at  150  m  is  needed  beyond  22  ms,  then  additional  observations 
are  necessary  at  distances  closer  than  110  m  to  the  explosion.  The  pressure 
history.  Figure  10,  shows  close  agreement  between  fitted  values  and  theoretical 
values.  The  magnitude  of  the  estimated  variances  of  the  pressure  is  of  the 
same  order  (5  percent)  as  the  assumed  accuracy  of  pressure  observations,  indi¬ 
cating  that  excessive  error  amplification  does  not  take  place. 


The  particle  velocity  history  is  shown  in  Figure  11.  The  figure  also 
contains,  as  a  dashed  line,  the  results  of  the  control  calculation  by  the 
quadrature  (2.15).  The  difference  between  both  velocity  curves  is  of  the  same 
order  as  the  estimated  standard  error  of  the  velocity.  Therefore,  one  can  con¬ 
clude  from  the  test  that  the  pressure  field  model  is  consistent  with  the  assumed 
accuracy  of  the  data.  Also  shown  in  the  figure  is  the  theoretical  particle 
velocity  as  predicted  by  the  strong  blast  theory.  The  deviation  between  the 
calculated  velocity  and  the  theoretical  velocity  is  always  in  the  same  direc¬ 
tion,  indicating  a  systematic  cause  for  the  difference.  A  systematic  reason 
is  indeed  present,  because  for  the  strong  blast  theory  one  sets  in  all  formulas 
the  ratio  pQ/ps  equal  to  zero,  whereas,  the  present  calculation  takes  the  ratio 
properly  into  account.  The  difference  in  the  result  can  be  easily  checked  at 
the  initial  point  of  the  velocity  history.  One  obtains  the  value  of  the  veloc¬ 
ity  as  predicted  by  strong  blast  theory  by  using  Eqs.  (4.3)  and  (4.4)  with 
Po/Ps  =  whereas  the  full  formulas  produce  the  value  shown  as  ’’calculated 
result.”  It  is  apparent  from  this  comparison  that  a  shock  pressure  ratio  of 
about  0.01  cannot  be  neglected  for  the  calculation  of  details  of  the  blast 
field . 


Figure  12  shows  the  density  history  at  150  m.  The  calculation  of  the 
density  makes  use  of  both  the  pressure  and  velocity  calculation  results. 
Nevertheless,  the  estimated  variance  of  the  result  is  moderate,  indicating  that 
the  computing  method  is  numerically  stable.  A  systematic  deviation  between  the 
calculated  density  and  the  strong  blast  density  is  again  present  and  caused  by 
the  neglect  of  the  ratio  pQ/ps  in  the  strong  blast  theory.  The  effect  of  the 
ratio  on  Eq.  (4.15)  is  stronger  for  lower  shock  pressures  and,  therefore,  the 
difference  is  larger  at  the  beginning  of  the  density  history,  where  the  parti¬ 
cles  have  been  subjected  to  a  less  strong  shock. 

The  dynamic  pressure  history  at  150  m  is  shown  in  Figure  13.  The  dynamic 
pressure  is  computed  from  the  velocity  and  density  and,  therefore,  one  obtains 
the  same  systematic  deviation  as  discussed  above  between  the  calculated  result 
and  the  strong  blast  theory.  More  interesting  is  the  estimated  variance  of 
the  calculated  dynamic  pressure.  It  is  of  the  order  of  five  percent,  i.e.,  of 
the  same  order  as  the  assumed  accuracy  of  the  pressure  data.  This  result  is 
very  encouraging  for  the  application  of  the  method  to  field  experiments,  be¬ 
cause  it  shows  that  the  computation  method  is  sufficiently  stable  to  produce 
accurate  dynamic  pressure  histories  from  overpressure  observations  with  typi¬ 
cal  uncertainties. 

Figures  14  through  17  show  the  calculated  field  histories  at  160  m  distance 
from  the  explosion.  The  results  are  similar  to  the  corresponding  results  for 
150  m  distance,  and  the  same  comments  apply  as  before  to  the  systematic  differ¬ 
ences  and  estimated  variances  of  the  calculated  histories. 
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PRESSURE 


BLRST  FIELD  BY  SBL-R0U TINES 
DISTRNCE  FROM  THE  EXPLOSION  150.00 
ERROR  LIMITS  CORRESPOND  TO  2 .00  STRNDRRD  ERRORS 


TIME 


Figure  10.  Pressure  History  at  150  m  in  Theoretical  Test  Example 


Field  fitting  results  with  confidence  limits 
Strong  blast  theory 


Time  is  expresses  in  seconds  and  pressure  is  expressed  in  pascals. 
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VELOCITY 


BLRST  FIELD  BY  SBL-ROUTINES 
DISTRNCE  FROM  THE  EXPLOSION  150.00 
ERROR  LIMITS  CORRESPOND  TO  2.00  STRNDQRD  ERRORS 


Figure  11.  Particle  Velocity  History  at  150  m  in  Theoretical 
Test  Example 


Calculated  result  with  confidence  limits 
Control  calculation 
Strong  blast  theory 


Time  is  expressed  in  seconds  and  velocity  is  expressed  in  m/s. 
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DENSITY 


BLAST  FIELD  BY  SBL-R0UTINES 
DISTANCE  FROM  THE  EXPLOSION  150.00 
ERROR  LIMITS  CORRESPOND  TO  2.00  STANDARD  ERRORS 


190  194  198  202  206  210  214  218x10^ 

TIME 


Figure  12.  Density  History  at  150  m  in  Theoretical  Test  Example 


Calculated  result  with  confidence  limits 
Strong  blast  theory 


Time  is  expressed  in  seconds  and  density  is  expressed  in  kg/m^. 
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DYNAMIC  PRESSURE  RH0*V**2/2 


BLAST  FIELD  BY  SBL-RQUTINES 
DISTANCE  FROM  THE  EXPLOSION  150.00 
ERROR  LIMITS  CORRESPOND  TO  2-00  STANDARD  ERRORS 


190  194  198  202  206  210  £14  218*10^ 

TIME 


Figure  13.  Dynamic  Pressure  History  at  150  m  in  Theoretical  Test  Example 


Calculated  result  with  confidence  limits 


I  j  Strong  blast  theory 


Time  is  expressed  in  seconds  and  dynamic  pressure  is  expressed  in  pascals. 


PRESSURE 


BLRST  FIELD  BY  SBL -ROUTINES 
DISTANCE  FROM  THE  EXPLOSION  160,00 
ERROR  LIMITS  CORRESPOND  TO  2-00  STANDARD  ERRORS 


22  23  24  25  26  27  *lOr* 

TIME 


Figure  14.  Overpressure  History  at  160  m  in  Theoretical  Test  Example 


Field  fitting  with  confidence  limits 
Strong  blast  theory 


Time  is  expressed  in  seconds  and  pressure  is  expressed  in  pascals. 
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VELOCITY 


BLRST  FIELD  BY  SBL -ROUTINES 
DISTANCE  FROM  THE  EXPLOSION  160.00 
ERROR  LIMITS  CORRESPOND  TO  2.00  STHNDRRD  ERRORS 


22  23  24  25  26  27  *10^ 

TIME 

Figure  15.  Particle  Velocity  History  at  160  m  in  Theoretical  Test  Example 


Calculated  result  with  confidence  limits 
Control  calculation 
Strong  blast  theory 


Time  is  expressed  in  seconds  and  velocity  is  expressed  in  m/s. 
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DENSITY 


BLRST  FIELD  BY  SBL -ROUTINES 
DISTRNCE  FROM  THE  EXPLOSION  160.00 
ERROR  LIMITS  CORRESPOND  TO  2.00  STRNDRRD  ERRORS 


22  23  2H  25  26  27  *10”3 

TIME 

Figure  16.  Density  History  at  160  m  in  Theoretical  Test  Example 


-4+h+H- 


Calculated  result  with  confidence  limit 
Strong  blast  theory 


Time  is  expressed  in  seconds  and 


density  is  expressed  in  kg/m 
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DYNAMIC  PRESSURE  RH0wV**2/2 


BLAST  FIELD  BY  SBL -ROUTINES 
DISTANCE  FROM  THE  EXPLOSION  160.00 
ERROR  LIMITS  CORRESPOND  TO  2.00  STANDARD  ERRORS 


22  23  2*4  25  26  27  *10^ 

TIME 

Figure  17 .  Dynamic  Pressure  History  at  160  m  in  Theoretical 
Test  Example 


-W-H-+4 


Calculated  result  with  confidence  limit 
Strong  blast  theory 


Time  is  expressed  in  seconds  and  dynamic  pressure  is  expressed  in  pascals. 
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8.  APPLICATION  TO  DATA  FROM  MISERS  BLUFF  II 

In  order  to  test  the  applicability  of  the  flow  field  calculation  method 
and  corresponding  computer  programs  to  real  problems,  we  used  the  computer  pro¬ 
grams  on  data  from  the  event  MISERS  BLUFF  11^. 

For  the  shock  fitting,  we  used  data  consisting  of  18  shock  arrival  times 
and  overpressure  measurements  at  ranges  between  89.7  m  and  715  m.  The  shock 
data  and  data  describing  the  ambient  conditions  are  listed  in  Table  5.  Table  6 
gives  the  numerical  results  of  the  shock  fitting.  Comparing  the  results  with 
corresponding  results  in  Table  3  for  the  theoretical  test  example,  one  notices 
that  now  there  are  no  indications  of  a  redundancy  of  parameters.  Their  esti¬ 
mated  standard  errors  are  relatively  small  and  correlations  between  them  are 
less  strong  than  in  the  theoretical  example.  The  standard  error  of  weight  one 
has  a  value  close  to  one,  indicating  that  the  estimated  standard  errors  of  the 
observations  are  consistent  with  the  shock  model  and  with  the  scatter  of  the 
data. 


Figures  18,  19,  and  20  show  the  shock  fitting  results  graphically.  The 
increase  of  the  estimated  confidence  interval  outside  the  fitting  region  is 
much  more  gradual  than  in  the  theoretical  test  example  (Figures  4,  5,  and  6). 
This  behavior  of  the  confidence  limits  can  be  attributed  to  the  fact  that  the 
model  contains  no  redundant  parameters  for  the  treatment  of  the  given  data  set. 

In  order  to  test  and  illustrate  the  flow  field  calculation  capability  of 
the  computer  program  we  assumed  that  the  stations  of  interest  were  located  at 
120  m  and  150  m  from  the  explosion;  i.e.,  we  assumed  that  the  goal  is  to  com¬ 
pute  flow  histories  at  these  distances.  Then  the  first  step  is  to  establish 
pressure  field  functions  in  the  vicinity  of  the  stations  of  interest.  Available 
were  overpressure  history  data  at  eight  stations  between  89.7  m  and  211  m,  as 
shown  in  Figure  21. 

The  overpressure  history  data  and  fitting  curves  for  the  individual  his¬ 
tories  are  shown  in  Figures  22  through  29-  For  this  data  fitting,  as  well  as 
for  the  subsequent  pressure  field  fitting,  we  assumed  the  following  standard 
errors  of  the  data: 

Pressure:  5%  of  the  initial  overpressure 

Time:  1  ms 

Range  and  Elevation:  1  m 


p 

G.  Teel,  "Free  Field  Airblast  Definition,  "  Proceedings  of  the  MISERS  BLUFF 
Phase  II  Results  Symposium ,  27-29  March  1979,  Defence  Nuclear  Agency  Report 
POR  7013-1 ,  1979. 
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TABLE  5 .  AMBIENT  AND  SHOCK  DATA  OF  MISERS  BLUFF  II 
Ambient  Conditions 

Molar  mass  M  =  28,94  g/mole 

Specific  heat  ratio  y  =  1.4 


Temperature 

T  =  312 

.05  K 

Pressure 

p  =  92. 
o 

92  kPa 

Charge 

Released  Energy  418.4  GJ 

Elevation 

0.0  m 

Shock 

Observations 

Range 

(m) 

Elevation 

(m) 

Overpressure 

(kPa) 

Time 

(ms) 

89.7 

0.0 

478 

56.6+1 

96.0 

0.0 

340 

64.8+1 

103.0 

0.0 

294 

75.0+1 

114 

0.0 

198 

92.4+2 

114 

0.0 

204 

95.7+2 

120.3 

0.0 

149 

104+2 

150.3 

0.0 

83.1 

164+2 

174 

0.0 

57.3 

217+2 

190 

0.0 

49.2 

254+2 

206 

1.22 

42.1 

292+2 

206 

2.44 

41.1 

293+2 

211 

0.0 

38.9 

304+3 

264 

0.0 

25.8 

436+3 

309 

1.22 

19.8 

551+3 

344 

0.0 

19.4 

643+3 

388 

0.0 

14.6 

760+4 

532 

0.0 

9.78 

1150+5 

715 

0.0 

6.59 

1660+10 

The  range  and  elevation  accuracies  were  assumed  to  be  1.0  m. 
The  overpressure  accuracy  was  assumed  to  be  5  percent. 
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TABLE  6.  SHOCK  PARAMETERS  OF  MISERS  BLUFF  II 


The  model  is  given  by  Eqs.  (4.7)  and  (4.8). 
The  data  are  listed  in  Table  5. 


Parameters  and  Standard  Errors 

a  =  (6.093  +  0.435) -106  Pa-m 
b  =  (-1.134  +  0. 217) *109  Pa-m2 
c  =  (3.343  +  0.235) -1011  Pa-m3 
d  =  57.934  +  0.766  ms  =  arrival  time  at  89.7  m 

The  standard  error  of  weight  one,  m0  =  1.156, 
is  included  in  the  quoted  standard  errors. 


Correlation  Matrix  of  the  Parameters 


(1.0 
-0.93948052 
0.83556376 
0.01536038 


-0.93948052 

1.0 

-0.96070389 

-0.03520487 


0.83556376 

-0.96070389 

1.0 

0.09994155 


0.01536038 

-0.03520487 

0.09994155 

1.0 


NOTE:  In  this  example,  the  standard  error  of  weight  one,  mQ,  is  close  to  one, 

indicating  that  the  data  scatter  is  consistent  with  the  data  standard  error 
as  specified  by  input.  However,  because  mc  is  larger  than  one,  the  inaccuracies 
might  have  been  slightly  underestimated.  We  have  multiplied  the  calculated 
standard  errors  of  the  parameters  by  m  in  order  to  account  for  the  possible 
underestimate.  ° 
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PRESSURE 


1 - 1 — T  1  1  TT 


DISTANCE 


MISERS  BLUFF  II 


CONFIDENCE  LIMITS  FOR  3.0  STANDARD  ERRORS 
WITHOUT  THE  FACTOR  ERZ  =  1.157 
ADJUSTED  ARE  OBSERVATIONS  OF 
PRESSURE,  DISTANCE  AND  TIME 


Figure  18.  Shock  Overpressure  Versus  Distance  in  MISERS  BLUFF  II 


Distance  is  expressed  in  meters  and  pressure  is  expressed  in  pascals. 
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MISERS  BLUFF  II 

CONFIDENCE  LIMITS  FOR  3.0  STANDARD  ERRORS 
WITHOUT  THE  FACTOR  ERZ  =  1 .157 
ADJUSTED  ARE  OBSERVATIONS  OF 
PRESSURE, DISTANCE  AND  TIME 


Figure  19.  Shock  Overpressure  Versus  Time  in  MISERS  BLUFF  IT 


Time  is  expressed  in  seconds  and  pressure  is  expressed  in  pascal 


DISTRNCE 


103 


102 


101 


TIME 


MISERS  BLUFF  II 

CONFIDENCE  LIMITS  FOR  3.0  STANDARD  ERRORS 
WITHOUT  THE  FACTOR  ER2  =  1.157 
ADJUSTED  ARE  OBSERVATIONS  OF 
PRESSURE, DISTANCE  AND  TIME 


Figure  20.  Shock  Distance  Versus  Time  in  MISERS  BLUFF  II 


Time  is  expressed  in  seconds  and  distance  is  expressed  in  meters. 
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Figure  21.  Locations  of  History  Observations  in  MISERS  BLUFF  II 


Distances  are  expressed  in  meters  and  times  are  expressed  in  seconds. 
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PRESSURE 


MISERS  BLUFF  II 

89.7  PS 

FITTED  CURVE  WITH  2-0  STRUGRRO  ERRORS 

*10  5 
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Figure  22.  Fitting  of  Pressure  at  89.7  m  in  MISERS  BLUFF  II 

The  data  accuracy  is  indicated  by  a  two-standard  error 
ellipse  around  each  observed  point. 

Time  is  expressed  in  seconds  and  pressure  is  expressed  in  pascals. 
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PRESSURE 


MISERS  BLUFF  II 
103  PS 

FITTED  CURVE  WITH  2-0  STfiNDP.R0  ERRORS 
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Figure  23.  Fitting  of  Pressure  at  103  m  in  MISERS  BLUFF  II 


The  data  accuracy  is  indicated  by  a  two-standard  error 
ellipse  around  each  observed  point. 


Time  is  expressed  in  seconds  and  pressure  is  expressed  in  pascals. 
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PRESSURE 
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FITTED  CURVE  WITH  P-0  STRNJCRO  ERRORS 
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Figure  24,  Fitting  of  Pressure  at  114  m  in  MISERS  BLUFF  II 


The  data  accuracy  is  indicated  by  a  two-standard  error 
ellipse  around  each  observed  point. 


Time  is  expressed  in  seconds  and  pressure  is  expressed  in  pascals. 
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PRESSURE 


MISERS  BLUFF  II 

IPO. 3  PS 

FITTED  CURVE  WITH  2-0  STFNDRRD  ERRORS 


TIME 


Figure  25.  Fitting  of  Pressure  at  120.3  m  in  MISERS  BLUFF  II 


The  data  accuracy  is  indicated  by  a  two-standard  error 
ellipse  around  each  observed  point. 


Time  is  expressed  in  seconds  and  pressure  is  expressed  in  pascals. 
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MISERS  BUJKF  [[ 
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KITTED  CURVE  WITH  2-0  STANDARD  ERRORS 
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Figure  26,  Fitting  of  Pressure  at  150.3  m  in  MISERS  BLUFF  XI 

The  data  accuracy  is  indicated  by  a  two-standard  error 
ellipse  around  each  observed  point. 

Time  is  expressed  in  seconds  and  pressure  is  expressed  in  pascals. 
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PRES5URE 


MISERS  BLUFF  II 
174  PS 


FITTED  CURVE  WITH  2-0  STRNDHRD  ERRORS 


TIME 

Figure  27.  Fitting  of  Pressure  at  174  m  in  MISERS  BLUFF  II 


The  data  accuracy  is  indicated  by  a  two-standard  error 
ellipse  around  each  observed  point. 


Time  is  expressed  in  seconds  and  pressure  is  expressed  in  pascals. 
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Figure  28.  Fitting  of  Pressure  at  206  m  in  MISERS  BLUFF  II 


The  data  accuracy  is  indicated  by  a  two-standard  error 
ellipse  around  each  observed  point. 


Time  is  expressed  in  seconds  and 


pressure  is  expressed  in  pascals. 
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Figure  29.  Fitting  of  Pressure  at  211  m  in  MISERS  BLUFF  II 


The  data  accuracy  is  indicated  by  a  two-standard  error 
ellipse  around  each  observed  point. 

Time  is  expressed  in  seconds  and  pressure  is  expressed  in  pascals. 


63 


Comparing  the  overpressure  histories,  one  notices  that  the  last  two  histories, 
Figures  28  and  29,  have  shapes  that  are  different  from  the  other  curves.  This 
observation  is  confirmed  by  the  plots  of  the  individual  history  fitting  param¬ 
eters  versus  distance,  Figure  30.  The  parameters  corresponding  to  the  last 
two  histories  do  not  follow  the  trends  of  the  other  parameters.  Based  on  this 
observation,  one  should  check  the  experimental  setup  and  see  if  any  causes  for 
the  different  behavior  of  the  last  two  histories  can  be  detected  that  may  ren¬ 
der  the  observations  less  reliable.  For  the  present  example,  we  assumed  that 
all  data  are  valid  and  equally  accurate. 

The  next  step  in  computing  the  flow  history  at,  say,  120  m  is  to  choose  a 
range  of  observations  for  the  overpressure  field  fitting.  The  range  should  be 
sufficiently  large  to  enable  one  to  calculate  the  flow  and  the  test  velocity 
for  the  time  interval  of  interest.  On  the  other  hand,  the  range  should  not  be 
too  large,  because  systematic  differences  between  the  model  and  the  solution 
of  the  governing  equations  may  become  excessive  in  a  large  fitting  region.  We 
illustrate  the  choice  of  a  fitting  region  by  using  two  regions  for  the  120  m 
history.  One  region  includes  all  histories  between  89.7  m  and  174  m,  and  the 
other  region  includes  histories  between  103  m  and  174  m.  (Figure  21  shows  that 
the  174  m  history  data  must  be  included,  at  least  partially,  in  order  to  com¬ 
pute  the  test  velocity  for  the  total  duration  of  the  120  m  observation.) 

Results  based  on  field  fitting  in  the  first  region  (89.7  m  -  174  m)  are 
shown  in  Figures  31  through  34.  The  dynamic  pressure,  which  is  an  important 
quantity  for  applications,  has  in  Figure  34  very  reasonable  estimated  confi¬ 
dence  limits.  However,  the  test  velocity  in  Figure  32  deviates  from  the  particle 
path  velocity  by  two  standard  errors  throughout  the  history.  Because  the  test 
velocity  is  computed  using  data  at  distances  beyond  120  m,  the  difference  in¬ 
dicates  a  discrepancy  between  data  below  120  m  and  beyond  120  m,  respectively. 

This  discrepancy  indication  disappears  if  one  uses  field  data  only  between 
103  m  and  174  m;  i.e.,  if  the  history  at  89.7  m  is  not  used.  The  results  of 
such  a  computation  are  shown  in  Figures  35  through  38.  The  velocity  plot  in 
Figure  36  now  shows  good  agreement  between  the  original  and  the  test  velocities. 
Comparing  Figure  36  with  Figure  32,  one  sees  that  the  test  velocity  in  both 
figures  is  practically  equal,  whereas  the  original  velocity,  which  depends  more 
directly  on  data  below  120  m,  has  changed  by  the  exclusion  of  the  89.7  m  his¬ 
tory.  Short  of  other  information  about  the  experiment,  one  should  in  such  a 
case  use  only  field  data  between  103  m  and  174  m  to  calculate  histories  at 
120  m,  because  then  the  results  are  self  consistent.  We  notice,  in  passing, 
that  the  absolute  differences  between  the  dynamic  pressures  in  both  calculated 
cases  are  small,  although  they  are  of  the  order  of  two  standard  errors. 

For  the  computation  of  the  flow  histories  at  150  m  range,  again  two  dif¬ 
ferent  pressure  field  functions  were  used.  One  function  was  obtained  by  fitting 
overpressure  data  from  ranges  between  120.3  m  and  206  m.  This  corresponds  to 
the  minimum  coverage  of  the  pressure  field  that  is  necessary  for  the  computation 
of  the  histories  at  105  m  for  the  complete  duration  of  the  pressure  observations. 
The  computed  histories  are  shown  in  Figures  39  through  42. 

Another  pressure  function  was  obtained  by  adding  two  more  observed  over¬ 
pressure  histories  to  the  data  base,  one  at  each  end  of  the  range  covered. 

Now,  the  overpressure  field  function  was  fitted  to  data  between  114  m  and  211  m. 
The  corresponding  histories  at  150  m  were  found  to  be  identical  within  the 
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PRESSURE 


MISERS  BLUFF  II/  (89.7-174.  M) 

DISTANCE  FROM  THE  EXPLOSION  120-00 
ERROR  LIMITS  CORRESPOND  TO  2.00  STANDARD  ERRORS 


TIME 


Figure  31  .  Pressure  History  at  120  m  in  MISERS  BLUFF  11/(89.7-174  m) 


Time  is  expressed  in  seconds  and  pressure  is  expressed  in  pascals. 
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VELOCITY 


MISERS  BLUFF  II/  (89.7-174.  M) 

DISTANCE  FROM  THE  EXPLOSION  120.00 
ERROR  LIMITS  CORRESPOND  TO  2.00  STANDARD  ERRORS 


Figure  32  .  Velocity  History  at  120  m  in  MISERS  BLUFF  11/(89.7-174  m) 


Velocity  and  confidenc 
Control  calculation 


e 


limits 


computed  along  particle  paths 


Time  is  expressed  in  seconds  and  velocity  is  expressed  in  m/s. 
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DENSITY 


MISERS  BLUFF  II/  C 89 -7- 174.  M) 

DISTRNCE  FROM  THE  EXPLOSION  120-00 
ERROR  LIMITS  CORRESPOND  TO  2-00  STANDARD  ERRORS 


Figure  33.  Density  History  at  120  m  in  MISERS  BLUFF  11/(89.7-174  m) 


Time  is  expressed  in  seconds  and  density  is  expressed  in  kg/m^. 
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DYNAMIC  PRESSURE  RH0*V**2/2 


MISERS  BLUFF  II/  (89.7-174.  M) 

DISTANCE  FROM  THE  EXPLOSION  120.00 
ERROR  LIMITS  CORRESPOND  TO  2.00  STANDARD  ERRORS 


9  10  11  12  13  14  15  16  17  18  19  20 

TIME 


Figure  34-.  Dynamic  Pressure  History  at  120  m  in  MISERS  BLUFF  11/(89.7-174  m) 


Time  is  expressed  in  seconds  and  dynamic  pressure  is  expressed  in  pascals. 
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PRESSURE 


.MISERS  BLUFF  II/  (103. -174.  M) 

DISTANCE  FROM  THE  EXPLOSION  120.00 
ERROR  LIMITS  CORRESPOND  TO  2.00  STANDARD  ERRORS 
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Figure  35.  Pressure  History  at  120  m  in  MISERS  BLUFF  11/(103-174  m) 


Time  is  expressed  in  seconds  and  pressure  is  expressed  in  pascals. 
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VELOCITY 


MISERS  BLUFF  II/  (103. -174.  M) 

DISTANCE  FROM  THE  EXPLOSION  120.00 
ERROR  LIMITS  CORRESPOND  TO  2.00  STANDARD  ERR  CHS 


Figure  36.  Velocity  History  at  120  m  in  MISERS  BLUFF  II/ (103-174  m) 


Velocity  and  confidence  limits  computed  along  particle  path 


^  Control  calculation 


Time  is  expressed  in  seconds  and  velocity  is  expressed  in  m/s. 
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DENSITY 


MISERS  BLUFF  II/  C 103.-174.  M) 

DISTANCE  FROM  THE  EXPLOSION  120.00 
ERROR  LIMITS  CORRESPOND  TO  2-00  STANDARD  ERRORS 
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Figure  37’.  Density  History  at  120  m  in  MISERS  BLUFF  11/(103-174  m) 


Time  is  expressed  in  seconds  and  density  is  expressed  in  kg/m  . 
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DYNAMIC  PRESSURE  RH0*V**2/2 


MISERS  BLUFF  11/  (103.-174.  Mi 
DISTANCE  FROM  THE  EXPLOSION  120.00 
ERROR  LIMITS  CORRESPOND  TO  2.00  STANDARD  ERRORS 


Figure  3S  .  Dynamic  Pressure  History  at  120  m  in  MISERS  BLUFF  11/(103-174  m) 
Time  is  expressed  in  seconds  and  dynamic  pressure  is  expressed  in  pascals. 
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PRESSURE 


MISERS  BLUFF  II/  C 120-3-206.  M 
DISTRMCE  FROM  THE  EXPLOSION  150-00 
ERROR  LIMITS  CORRESPOND  TO  2-00  STRNDRRD  ERRORS 


15  17  19  21  23  25  27  aTiD--2 

TIME 

Figure  39.  Pressure  History  at  150  m  in  MISERS  BLUFF  11/(120.3-206  m) 


Time  is  expressed  in  seconds  and  pressure  is  expressed  in  pascals. 
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VELOCITY 


MISERS  BLUFF  II/  (120.3-206.  M 
DISTANCE  FROM  THE  EXPLOSION  L50.00 
ERROR  LIMITS  CORRESPOND  TO  2.00  STRNDfiRD  ERRORS 
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TIME 

Figure  40.  Velocity  History  at  150  m  in  MISERS  BLUFF  11/(120.3-206  m) 


Velocity  and  confidence  limits  computed  along  particle  paths 


Control  calculation 


Time  is  expressed  in  seconds  and  velocity  is  expressed  in  m/s. 
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DENSITY 


MISERS  BLUFF  II/  (120.3-206.  M 
DISTANCE  FROM  THE  EXPLOSION  150.00 
ERROR  LIMITS  CORRESPOND  TO  2.00  STANDRRO  ERRORS 


15  17  19  21  23  25  27  *10^ 

TIME 

Figure  41.  Density  History  at  150  m  in  MISERS  BLUFF  11/(120.3-206  m) 

O 

Time  is  expressed  in  seconds  and  density  is  expressed  in  kg/in  . 
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DYNAMIC  PRESSURE  RH0*V**2/2 


MISERS  BLUFF  II/  (120.3-206.  M 
DISTANCE  FROM  THE  EXPLOSION  150.00 
ERROR  LIMITS  CORRESPOND  TO  2.00  STANDARD  ERRORS 


15  11  L9  21  23  25  27  xtO-2 

TIME 


Figure  42.  Dynamic  Pressure  History  at  150  m  in  MISERS  BLUFF  11/(120.3-206  m) 


Time  is  expressed  in  seconds  and  dynamic  pressure  is  expressed  in  pascals. 
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plotting  accuracy  to  those  of  the  previous  calculation.  A  closer  inspection 
of  the  numerical  output  indicates  that  the  variance  estimates  are  slightly 
reduced  by  the  addition  of  the  two  observed  histories.  However,  the  difference 
between  streamline  velocity  and  test  velocity  is  increased.  The  latter  obser¬ 
vation  confirms  the  experience  with  history  calculations  at  120  m,  where  it  was 
also  found  that  best  results  are  obtained  when  the  range  for  the  fitting  of  the 
pressure  field  function  is  as  small  as  practically  possible. 

Of  some  concern  in  the  present  example  is  the  large  difference  between 
the  original  and  test  velocity,  which  is  of  the  order  of  two  standard  errors. 

This  difference  may  be  attributed  to  disturbances  in  the  flow  field  at  ranges 
beyond  150  m,  that  can  violate  some  of  the  basic  assumptions  for  the  flow  model, 
e.g.,  the  spherical  symmetry.  In  a  practical  application,  one  would  reexamine 
the  data  and  then  decide,  depending  on  the  purpose  of  the  history  calculations, 
whether  the  proposed  analysis  is  adequate. 

The  example  makes  it  clear  that  one  needs  a  good  definition  of  the  pres¬ 
sure  field  around  the  reference  ground  range  where  the  flow  history  is  of 
interest.  The  range  interval  in  which  additional  pressure  histories  should  be 
recorded  depends  on  the  size  of  the  explosion,  on  the  reference  range,  and  on 
the  anticipated  duration  of  the  history  observation.  For  a  planned  experiment, 
one  may  assume  that  typically  the  duration  of  interest  will  approximate  the 
duration  of  positive  overpressure.  Estimates  of  positive  overpressure  duration 
are  available  from  past  experiments.  Figure  43  shows  such  an  estimate^  for 
an  0.1  kton  TNT  experiment,  like  MISERS  BLUFF. 

Estimates  of  corresponding  range  intervals  can  be  calculated,  e.g.,  assuming 
that  the  shock  trajectory  and  the  particle  paths  are  straight  lines  within  the 
area  of  interest.  Range  intervals  computed  with  this  assumption  are  overesti¬ 
mates,  as  can  be  seen  from  Figure  21.  The  results  are  as  follows. 

The  downstream  range  interval  that  is  needed  to  carry  out  the  control  cal¬ 
culation  for  a  history  of  duration  At  can  be  estimated  by 


Ar 


control 


At  •  U 


(8.1) 


where  U  is  the  shock  velocity.  The  corresponding  upstream  interval  that  is 
needed  to  compute  the  flow  field  history  at  the  reference  range  for  the  duration 
At  can  be  estimated  by 


Ar 

computing 


=  At (1/ u 

s 


-  1/U) 


At  *  U 


(8.2) 


where  ug  is  the  particle  velocity  behind  the  shock.  Using  the  shock  formula 
(4.4)  one  can  also  express  Eq.  (8.2)  in  the  form 


1  Ps 

Ar  .  =  At  •  U  -  — 

computing  y  p 


2y  PQ 


(8.3) 


10 


Klaus  Opalka,  Terminal  Ballistics  Division,  USA  AREADCOM/Ballistic  Research 
Laboratory,  private  communication,  March  1981. 
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TIME  INTERVAL  (ms) 


Figure  43.  Estimated  Positive  Overpressure  Duration  in  a 
0.1  kton  TNT  Experiment 
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where  pg  is  the  shock  overpressure  and  p0  is  the  ambient  pressure.  The  esti¬ 
mates  (8.1)  and  (8.3)  of  the  downstream  and  upstream  range  intervals  are  plotted 
in  Figure  44  as  functions  of  the  reference  range,  using  the  At  from  Figure  43. 

In  order  to  check  the  quality  of  the  approximate  range  interval  formulas, 
we  also  have  determined  the  range  intervals  by  an  analysis  of  Figure  21.  The 
resulting  upstream  intervals  are  plotted  in  Figure  44.  The  plot  shows  that 
for  close  ranges  the  approximate  formulas  overestimate  the  range  interval  by 
more  than  a  factor  of  two.  For  ranges  larger  than  200  m  the  factor  is  smaller 
because  the  path  line  curvatures  decrease  with  increasing  range. 

The  downstream  range  intervals  computed  by  an  analysis  of  Figure  21  differ 
only  insignificantly  from  the  curve  in  Figure  44  and,  therefore,  are  not  plotted. 


9.  CONCLUSIONS 

A  distinct  feature  of  the  described  method  for  blast  field  reconstruction 
from  overpressure  measurements  is  the  calculation  of  variance  estimates  of  all 
field  variables  concurrently  with  the  calculation  of  the  variables  themselves. 

The  estimates  can  be  used  in  the  planning  of  blast  experiments;  e.g.,  to  deter¬ 
mine  pressure  gage  arrangements  that  enhance  the  predicted  accuracy  of  the 
calculated  flow  field.  In  examples  where  the  method  was  applied  to  actual 
measurements,  it  was  found  that  the  sensitivity  of  the  computed  results  to 
observational  inaccuracies  was  up  to  ten  times  smaller  than  that  of  some  pre¬ 
vious  methods. 

Another  unique  feature  of  the  method  is  the  control  calculation  of  a  flow 
velocity  based  on  the  continuity  equation.  The  control  calculation  provides 
an  indicator  for  the  consistency  of  the  data  with  the  basic  assumption  for  the 
mathematical  model.  Numerical  experimentation  with  this  indicator  shows  that 
more  consistent  results  can  be  expected  when  the  pressure  field  function  is 
fitted  to  data  covering  smaller  ranges  of  distances.  This  means  for  the  planning 
of  experiments  that  one  should  provide  a  sufficient  number  of  pressure  history 
observations  in  the  vicinity  of  the  range  of  interest.  Favorable  arrangements 
of  stations  and  durations  of  pressure  recordings  can  be  obtained  from  test  runs 
of  the  field  calculation  programs  on  approximate  expected  data  sets. 

The  nine  parameter  pressure  field  function  of  this  report  was  found  to  be 
sufficiently  flexible  for  the  description  of  the  pressure  within  a  wide  range 
of  theoretical  and  real  problems. 

One  part  of  the  analysis,  having  independent  applications,  is  the  deter¬ 
mination  of  shock  functions  from  shock  data.  The  method  described  in  this 
report  has  the  advantage  of  using  all  available  data,  including  incomplete 
sets  and  adjusting  simultaneously  measurements  of  pressure,  time,  and  distance. 
The  method  can  be  expected  to  produce  more  reliable  results  than  previous  nu¬ 
merical  shock  fittings. 
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RANGE  INTERVALS  (m) 


Figure'  44.  Range  Intervals  for  a  0.1  kton  TNT  Experiment 
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LIST  OF  SYMBOLS 


(a,b,c,d)  -  shock  parameters 

(A,B,C)  -  overpressure  field  parameter  functions 
(A^AgjB^Bg,^)  -  overpressure  field  parameters 
c  -  sound  speed  (m/s) 
e  -  specific  internal  energy  (J/kg) 

E  -  function  defined  by  Eq.  (5.3) 

f  -  model  function  for  data  fitting 

H  -  flow  field  vector  defined  by  Eq.  (3.17) 

M  -  molar  mass  (kg/mole) 
p,P  -  pressure  (Pa) 

p^  -  overpressure  field  function  (Pa) 
p^  -  overpressure  history  function  (Pa) 

Q  -  function  defined  by  Eq.  (4.16) 

r,R  -  distance  from  the  center  of  explosion 

t,T  -  time  after  explosion  (s) 

Tq  -  ambient  temperature  (K) 
u  -  particle  velocity  (m/s) 

U  -  shock  velocity  (m/s) 

V  -  variance-covariance  matrix 
y  -  ratio  of  specific  heats 

T  =  (y+1)/(2ypq)  -  constant  in  Section  4  (Pa”^-) 

0  -  parameter  vector  for  model  fitting 

3 

p  -  density  (kg/m  ) 

SUBSCRIPTS 

s  -  pertaining  to  shock 
o  -  pertaining  to  ambient  conditions 
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Department  MIS  42-37 
ATTN:  R.M.  Schmidt 
P.0.  Box  3999 
Seattle,  WA  98124 

1  General  Electric  Company-TEMPO 

ATTN:  DAS  LAC 
P . 0 .  Drawer  QQ 
Santa  Barbara,  CA  93102 

1  General  Electric  Company-TEMPO 

ATTN:  E.  Bryant 
220  S.  Main  Street,  Rm  206 
Bel  Air,  MD  21014 

2  Kaman  Avidyne 
ATTN:  N.P.  Hobbs 

S.  Criscione 
83  Second  Avenue 
Northwest  Industrial  Park 
Burlington,  MA  01830 

2  Kaman  Sciences  Corporation 

ATTN:  Library 

F.H.  Shelton 
Suite  855 

101  Continental  Boulevard 
El  Segundo,  CA  90245 

1  Pacific  Sierra  Research  Corp. 

ATTN:  H.  Brode 

1456  Cloverfield  Boulevard 
Santa  Monica,  CA  90404 

1  Physics  International  Corp. 

ATTN:  F.M.  Sauer 

J.M.  Thomsen 
Technical  Library 
2700  Merced  Street 
San  Leandro,  CA  94577 
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No.  of 
Copies 


DISTRIBUTION  LIST 


Organization 


No.  of 

Copies  Organization 


1  R&D  Associates  2 

ATTN:  H.  Cooper 

Suite  500 

1401  Wilson  Boulevard 
Arlington,  VA  22209 

6  R&D  Associates  3 

ATTN :  J .  Carpenter 

A.L.  Kuhl 
A.L.  Latter 
J.G.  Lewis 
W.B.  Wright 
Technical  Library 

P.0.  Box  9695  2 

2  The  Rand  Corporation 

ATTN :  C . C .  Mow 

Library-D 
1700  Main  Street 

Santa  Monica,  CA  90406  1 

2  Science  Applications,  Inc. 

ATTN:  J.  Cockayne 

W'.  Layson 

1710  Goodridge  Drive 

P.0.  Box  1303  2 

McLean,  VA  22101 

1  Science  Applications,  Inc. 

2450  Washington  Avenue 
San  Leandro,  CA  94577 

1 

2  Science  Applications,  Inc. 

ATTN:  Technical  Library 

M .  McKay 
P.0.  Box  2351 
La  Jolla,  CA  92037 

2 

2  Science  Applications,  Inc. 

ATTN:  J.F.  Dishon,  III 
R.I.  Miller 
Albuquerque,  NM  87115 

2 


Systems,  Science  and  Software 
ATTN :  C . E .  Needham 

J.E.  Crepeau 
P.0.  Box  8243 
Albuquerque,  NM  87198 

Terra  Tek,  Inc. 

ATTN :  S .  Green 

A.H.  Jones 
Technical  Library 
420  Wakara  Way 
Salt  Lake  City,  UT  84108 

TRW  Systems  Group 
ATTN:  B.  Sussholts 

Tech  Info  Ctr/S-1930 
One  Space  Park 
Redonodo  Beach,  CA  92078 

Weidlinger  Associate 
Consulting  Engineers 
ATTN:  M.L.  Baron 

110  East  59th  Street 
New  York,  NY  10022 

Drexel  Institute  of  Technology 
ATTN :  P . C .  Chou 
Mechanics  and  Structures  Group 
32  and  Chestnuts  Streets 
Philadelphia,  PA  19137 

Massachusetts  Institute  of  Technology 
Aeroelastic  and  Structures 
Research  Laboratory 
ATTN:  E.A.  Witmer 
Cambridge,  MA  02139 

Southwest  Research  Institute 
ATTN:  W.E.  Baker 
A.B.  Wenzel 
8500  Culebra  Road 
San  Antonio,  TX  78206 

SRI  International 
ATTN:  G.R.  Abrahamson 

C.  Peterson 
333  Ravenswood  Avenue 
Menlow  Park,  CA  94025 
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No.  of 
Copies 


DISTRIBUTION  LIST 


Organization 

2  University  of  Denver 

Denver  Research  Institute 
ATTN :  J .  Wisotski 

Technical  Library 
P.0.  Box  10127 
Denver ,  CO  80210 

1  J.  D.  Haltiwanger 

Consulting  Services 
B106a  Civil  Engineering  Bldg 
208  N.  Roraine  Street 
Urbana,  IL  61801 


No.  of 

Copies  Organization 

1  University  of  Washington 

ATTN:  K.A.  Holsapple 
Seattle,  WA  98195 

Aberdeen  Proving  Ground 

Dir,  USAMSAA 
ATTN:  DRXSY-D 

DRSXY-MP ,  H.  Cohen 
Cdr ,  USATECOM 

ATTN:  DRSTE-TO-F 
Dir,  USACSL,  Bldg.  E3516,  EA 
ATTN:  DRDAR-CLB-PA 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name: 


Telephone  Number: 
Organization  Address: 


